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In  this  investigation  a  mathematical  model  of  a  continuous 
stage-wise  distillation  column  was  derived.   A  convergence  technique  was 
developed  to  solve  the  mathematical  model  with  the  use  of  a  digital 
computer.   Experimental  data  from  an  industrial  column  was  used  to 
verify  the  accuracy  of  the  model  to  simulate  plant  conditions. 

The  mathematical  model  was  then  modified  so  that  it  could  be 
used  to  predict  the  steady-state  control  points  to  maintain  composition 
quality  in  one  of  the  product  streams  of  this  column.   This  control 
model  was  then  used  to  control  a  tower  in  a  tall  oil  production  unit. 

The  mathematical  model  consists  of  steady-state  material  and 
energy  balances  around  the  upper  and  lower  sections  of  the  fractionating 
column  after  the  method  of  Sorel.   The  convergence  method  uses  a  ratio 


of  feed  plate  compositions,  as  computed  from  the  top  and  bottom,  raised 
to  a  power  in  order  to  adjust  the  top  and  bottom  compositions  for  the 
next  trial.   Convergence  is  reached  when  the  ratio  equals  one.   The 
control  model  convergence  technique  adjusts  overheads  flow  rate  in 
order  to  maintain  constant  composition  in  the  sidestream  during  the 
convergence  procedure. 

An  IBM  1620/1070  computer  system  was  programmed  with  the  control 
model  and  linked  to  an  industrial  tall  oil  fatty  acid  fractionating 
column.   Vapor  phase  chromatograph  analyses  of  the  feed  stream  were 
entered  into  the  computer  at  regular  intervals  and  the  model  then 
predicted  external  flow  rates  necessary  to  maintain  constant  product 
composition  in  the  sidestream.   A  1  per  cent  increase  in  yield  of 
the  main  product  resulted  from  this  control.   This  would  produce  a 
financial  return  of  $10,000  a  month  for  the  average  tall  oil  plant  with 
a  throughput  of  10,000  tons  per  month. 

On  the  basis  of  this  work,  several  conclusions  may  be  made: 
It  is  possible  to  write  a  theoretical  description  of  an  actual 
industrial  distillation  column  with  sufficient  accuracy  to  allow  its 
control  by  a  digital  computer. 

It  is  not  always  necessary  to  have  all  physical  data  completely 
defined  before  a  model  may  be  used,  but  actual  operating  data  may  be 
used  to  fill  the  gaps  in  our  knowledge. 

The  power  ratio  convergence  technique  ha's  proved  itself  to  be 
a  useful  method  in  the  solution  of  steady-state  distillation  calculations 
and  should  have  applications  in  other  types  of  trial  and  error  calcula- 
tions . 


CUAriER  I 

INfRODUCTION 
A.   The  Process  of  Dlstilldtion  and  Its  Amenability  to  Control 

Distillation  is  one  of  the  major  unit  operations  in  the  chemical 
and  related  industries.   Entire  plants  consist  of  distillation  towers 
and  their  auxiliary  equipment.   The  control  of  these  tov^'ers  and  plants 
has  occupied  many  working  hours  of  many  process  chemical  engineers. 

Before  World  War  II,  process  control  was  more  of  an  art  than 
a  science.   It  consisted  mainly  of  measurement  and  control  of  various 
parts  of  the  process  independently  of  each  other  and  with  little  regard 
for  the  overall  process  involved.   Control  was  regarded  as  primarily  a 
function  of  instruments  and  hardware  (I)  and  not  of  the  process.   With 
the  development  of  an  adequate  theory  of  control,  the  possibility  of 
integrating  the  process  and  its  controllers  into  one  overall  system  for 
analysis  and  design  became  a  practical  goal,  and  much  effort  in  recent 
years  has  been  in  this  direction.   This  approach  made  necessary  considera- 
tion of  the  process  dynamics  and  the  realization  of  the  need  of  rapid 
and  reliable  continuous  composition  analyzers.   Though  much  work  is 
still  to  be  done  in  this  area,  considerable  progress  is  evident  from  the 
results  being  reported  in  the  literature. 

One  of  the  first  things  one  realizes  when  trying  to  effect 
contro]  of  a  unit  operation  is  that  the  two  concepts  of  design  and 
control  have  in  the  past  been  thought  of  as  separate  entities,  whereas 
good  design  should  include  consideration  of  dynamics  and  control.   Classi- 
cal design  equations  have  limited  utility  in  predicting  how  a  plant  will 
operate  (1).   Design  is  always  concerned  with  size  and  operating  conditions 
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are  assumed  to  be  fixed.   But  once  the  design  is  implemented,  process 
people  want  and  have  to  know  how  the  plant  operates  under  a  variety  of 
conditions,  regardless  of  what  it  was  designed  to  do.   Design  equations 
and  procedures  for  the  sizing  of  equipment  do  not  lend  themselves 
readily  to  the  prediction  of  operating  conditions  without  the  use  of 
trial  and  error  procedures.   However,  it  is  reasonable  to  assume  that 
the  control  problem  can  always  be  formulated  in  a  way  which  will  yield 
answers  to  the  questions  of  how  the  plant  operates. 

The  study  of  an  operating  system  can  be  thought  of  in  two  ways, 
which  while  related,  may  be  treated  as  separate  topics.   The  process 
may  be  described  as  a  time- invariant  operation,  or  in  a  steady  state, 
or  a  more  complete  description  of  the  process  can  be  formulated  in  which 
time  varying  conditions  are  included.   The  latter  is  the  dynamic  model 
approach  and  of  course  includes  the  steady-state  information  as  a  special 
case. 

1.   Steady  State 

Most  design  is  carried  out  on  a  steady- state  basis.   That  is, 
the  equipment  is  sized  to  produce  certain  results  when  operating  at 
constant  conditions.   Moreover,  most  industrial  operations  are  geared  for 
operation  at  steady-state  conditions  because  the  process  engineers  and 
operators  want  to  "line  out"  their  column  and  hold  the  product  on  grade 
to  meet  sales  specifications.   Thus,  it  is  not  surprising  to  find  many 
applications  for  control  of  a  distillation  tower  on  a  steady-state  basis. 
In  addition,  the  added  time  and  expense  to  describe  the  system  dynamically 
may  be  difficult  to  justify  from  an  economic  point  of  view,  even  though 
it  would  appear  obvious  that  a  thorough  dynamic  analysis  may  yield 


dividends  in  processing  knowledge  that  would  not  be  obtainable  from  a 
steady  state  consideration. 

The  steady-state  problem  can  be  stated  in  this  fashion:  Given 
a  feed  rate  and  composition,  at  what  values  must  the  controlled  varia- 
bles be  set  so  as  to  achieve  "best  operation"  (i.e.,  on-grade  product) 
when  the  column  reaches  a  steady  state? 

2 .   Dynamic 


Dynamic  distillation  models  have  been  formulated  by  several 
investigators.   Huckaba ,  Franke,  and  May  (2)  described  a  model  for 
binary  distillation  and  Distefano  (3)  used  this  model  in  an  open  loop 
predictive  control  scheme  on  a  pilot  plant  size  laboratory  column  with 
good  success.   Other  dynamic  models  have  appeared  in  the  literature 
(4,5,6),  but  their  application  to  control  has  been  limited.   It  is 
reasonable  to  assume  that  some  operating  advantage  could  be  obtained  by 
using  transient  conditions  to  operate  a  better  column;  however,  no 
published  description  of  such  an  operation  has  come  to  the  attention  of 
the  author.   Some  systems,  such  as  batch  distillation,  are  inherently 
dynamic  and  the  dynamic  approach  is  necessary  to  describe  and  understand 
them. 

In  order  to  use  a  dynamic  model  in  an  industrial  control  problem, 
and  not  in  the  study  of  the  system  for  design  of  the  control,  several 
conditions  must  be  present.   Perhaps  the  first  one,  which  is  basic  to  all 
engineering  problems,  is  that  the  financial  return  for  the  more  complex 
dynamic  study  must  be  great  enough  to  pay  out  the  additional  expense. 
There  is  no  general  answer  to  this  problem  either  one  way  or  the  other, 
but  the  largest  drawback  to  a  dynamic  analysis  at  the  present  time  is  the 


lack  of  experience  in  these  sorts  of  techniques.   Another  consideration 
is  the  availability  of  physical  data  and  the  obtaining  of  experimental 
verifications  of  the  dynamic  model.   In  the  plant  project  discussed  in 
this  work,  the  difficulty  of  obtaining  samples  from  the  process  and  the 
long  involved  chemical  analysis  of  these  samples  made  the  verification 
of  any  dynamic  model  impractical.   A  third  consideration  is  the  estimated 
response  time  of  the  process  and  the  frequency  and  magnitude  of  process 
upsets.   If  the  process  has  a  very  fast  response  time,  then  the  dynamic 
considerations  are  perhaps  not  as  important  as  the  more  usual  slow 
response  that  is  found  in  the  chemical  process  industry.   Also,  if 
upsets  are  infrequent,  then  the  dynamic  problem  is  a  minor  consideration. 
If  the  upsets  are  low  in  magnitude  and  slowly  changing,  rather  than 
abruptly,  it  would  be  possible  to  handle  the  dynamics  by  alert  operators 
rather  than  a  computer.   One  last  consideration  is  the  available  computer 
power.   Dynamic  analysis  requires  a  high  speed  digital  computer  to  be 
used  in  an  on-line  basis.   In  addition,  even  an  off-line  study  requires 
large  core  and  very  fast  computational  speed  to  provide  answers  in  a 
reasonable  amount  of  time.   Again,  in  the  plant  project,  we  were 
restricted  to  the  use  of  an  IBM  1620  with  limited  access  to  the  University 
of  Florida's  IBM  709.   The  IBM  1620  is  a  medium  speed  computer  and  takes 
in  the  range  of  one  to  two  hours  to  solve  the  steady- state  analysis  of 
the  distillation  tower. 

B.   Previous  Work 
The  literature  of  chemical  engineering  is  replete  with  methods 
to  calculate  the  performance  of  distillation  towers.   As  pointed  out  in 
Part  A,  most  of  these  techniques  are  design  techni-iues  which  give  size 
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and  amount  of  equipment  rather  than  predict  the  performance  of  existing 
equipment . 

The  complexity  of  the  equations  and  the  large  number  of  variables 
in  any  distillation  calculation  make  the  job  of  calculating  the  conditions 
within  the  column  a  formidable  one.   However,  only  a  relatively  small 
number  of  these  variables  are  independent  and  thus  have  to  be  specified 
before  the  problem  is  determinate.   The  number  of  specifications  and 
the  way  these  should  be  picked  have  been  discussed  by  Gilliland  and 
Reed  (7),  Kwauk  (8),  and  others  (9,10).   An  analysis  of  the  degrees  of 
freedom  for  this  problem  is  discussed  in  Appendix  A. 

1.   Calculation  Procedures  for  the  Steady  State 

a.   Short-cut  Methods 


Short-cut  methods  which  give  approximate  solutions  to  the 
rigorous  distillation  equations  continue  to  serve  a  useful  purpose  even 
with  the  availability  of  high-speed  computers  which  reduce  the  work 
necessary  to  obtain  rigorous  solutions.   Often  the  data  available  do 
not  warrant  the  use  of  rigorous  equations,  or  perhaps  a  large  number  of 
different  design  configurations  should  be  examined  before  attention  is 
focused  on  the  detailed  analysis  of  one  column. 

Fenske  (11)  has  presented  an  equation  which  relates  the  separa- 
tion obtained  between  tv;o  components  at  total  reflux  to  the  number  of 
equilibrium  stages.   Although  the  equation  is  derived  for  a  binary,  it 
may  be  used  for  a  pair  of  components  in  a  mult  1- component  mixture. 
Underwood  (12,13,14,15)  developed  eciuations  which  are  convenient  for 
estimating  the  minimum  reflux  corresponding  to  any  specified  separation 
of  two  key  components  in  a  mult icomponent  feed. 


Gilliland  (16)  has  presented  a  useful  correlation  which  relates 
minimum  stages  (calculated  with  Fenske  equations)  and  minimum  reflux 
(calculated  with  Underwood  equations)  to  the  actual  number  of  equilibrium 
stages  and  the  actual  reflux  ratio. 

Sauer  (17)  has  presented  a  computer  program  which  uses  the 
Fenske  and  Underwood  equations  and  makes  use  of  the  Gilliland  correla- 
tions to  calculate  theoretical  stages  and  reflux  ratio. 

Smith  (9)  presents  a  general  short-cut  method  which  is  applicable 
to  all  equilibrium  stage  processes  and  he  recommends  its  use. 

It  should  be  pointed  out  that  there  are  many  simplifying  assump- 
tions with  any  of  the  short-cut  methods.   The  assumption  of  optimum 
feed  location  is  inherent  in  the  Underwood  equations;  the  assumption  of 
saturated  reflux  is  inherent  in  both  the  Fenske  and  Underwood  equations. 
Constant  molal  overflow  is  assumed  in  all  these  short-cut  methods.   To 
use  the  Fenske  equations,  it  is  necessary  to  estimate  an  average  relative 
volatility  for  each  component.   For  very  small  volatility  differences, 
the  Fenske  equations  can  be  widely  in  error  (9). 

Harbert  (18)  also  presented  equations  for  the  case  of  constant 
molal  overflow  and  constant  relative  volatilities. 

b.  Rigorous  Methods 

Sorel  (19)  developed  and  applied  the  mathematical  theory  of  the 
rectifying  column  for  binary  mixtures.   By  making  energy  and  material 
balances  around  each  plate  and  assuming  that  equilibrium  was  attained 
between  the  vapor  and  liquid  leaving  the  plate,  he  was  able  to  calculate 
stepwise  through  the  column.   All  of  the  more  modern  methods  have  been  bui 
upon  Sorel 's  foundation. 


Lewis  and  Matheson  (20)  applied  Sorel's  method  to  multicoitiponent 
mixtures.   As  originally  presented,  their  method  was  concerned  with 
picking  the  number  of  ideal  stages,  feed  location  and  operating  pressure 
for  the  column.   Robinson  and  Gilliland  (21)  present  a  discussion  of  this 
method,  and  Bonner  (22)  has  modified  the  Lewis  and  Matheson  method  and 
made  it  useful  for  machine  computation.   Smith  (9)  gives  a  detailed 
discussion  of  Bonner's  use  of  the  method. 

The  other  m.ost  common  technique  for  solving  multicomponent 
distillation  problems  is  the  Thiele-Geddes  (23)  method.   This  method 
differs  from  the  Lewis  and  Matheson  method  in  that  rather  than  assuming 
top  and  bottom  compositions,  a  complete  temperature  profile  for  the 
column  is  assumed.   These  two  methods  differ  only  in  computational 
techni4ue;  the  final  results  of  the  calculations  are  the  same.   Some 
authors  have  shovm  distinct  preference  for  one  technique  or  the  other. 
Holland  (24)  states  that  his  convergence  method  works  best  with  the 
Thiele-Geddes  method.   The  choice  of  which  method  should  be  used  depends 
to  a  large  extent  on  the  particular  problem  being  solved.   The  conver- 
gence method  presented  in  this  work  is  independent  of  the  distillation 
calculations  and  works  equally  well  with  both  Lewis-Matheson  and 
Thiele-Geddes  techniques. 

Other  multicomponent  distillation  calculation  procedures  include 
that  of  Amundson-Pontinen  (25,26),  in  which  the  equations  are  written  as 
a  system  of  simultaneous  algebraic  equations  and  solved  by  matrix 
techniques.   Rose,  et  al.  (27)  have  applied  a  relaxation  technique  to 
carry  the  solutions  of  the  differential-difference  equations  to  steady 
state.   This  is  slower  than  most  other  methods  for  obtaining  the  steady- 


state  solution,  but  has  the  advantage  of  allowing  one  to  follow  the 
transient  response  from  the  disturbance  to  steady  state  if  realistic 
initial  conditions  are  specified.   Ball  (28)  has  modified  the  Rose 
relaxation  method  and  speeded  it  up  to  make  it  competitive  with  other 
methods  for  obtaining  steady-state  solutions  at  the  expense  of  correct 
dynamic  results,   Rosenbrock  (29)  has  also  presented  a  method  for  obtain- 
ing rapid  steady-state  solutions  from  a  dynamic  model.   With  the  advent 
of  better  and  more  efficient  numerical  methods,  steady-state  solutions 
to  a  system  of  differential-difference  equations  may  well  prove  the  most 
practical  method  to  make  multicomponent  distillation  calculations. 
This  will  undoubtedly  be  tested  as  the  transient  response  becomes 
more  common  in  its  use  by  design  and  operating  personnel. 

2.   Control  Applications 

Many  articles  have  been  written  on  the  control  of  distillation 
towers;  however  few  take  the  view  that  it  is  necessary  to  have  a  rigorous 
mathematical  model  to  describe  the  process.   Most  attempts  at  control 
have  focused  on  a  particular  problem  (30)  and  the  development  of  hardware 
to  handle  a  particular  job.   Buckley  (1)  in  his  book,  Techniques  of 
Process  Control,  devotes  four  chapters  exclusively  to  control  of  distilla- 
tion towers.   The  reader  is  referred  to  this  work  as  a  scholarly 
exposition  of  the  state  of  the  art. 

Horn  and  Miller  (31)  have  described  the  fitting  of  "an  available 
complex,  plate-to-plate  engineering-calculation  computer  program"  to 
actual  plant  data.   They  use  this  model  to  generate  data  for  a  regression 
analysis  to  obtain  simple  approximate  relations  among  the  control 
variables . 


MacMullan  and  Shinsky  (32)  report  the  feedforward  control  of  a 
super  fractionator  using  a  pneumatic  analog  controller.   They  use  a 
material  balance  equation  over  the  entire  column  for  steady  state  and 
postulate  a  second-order  model  with  one  lead  and  tvro  lag  terms  to 
simulate  the  dynamics  of  the  system.   They  report  substantial  improvement 
over  the  more  conventional  feedback  control  techniques,  reducing  the 
distillation  composition  variation  from  5  per  cent  to  1  per  cent. 

C.   Statement  of  the  Problem 

The  problem  of  digital  computer  control  of  distillation  processes 
is  a  general  one  but  its  realistic  application  has  to  be  directed  at  a 
particular  process.   The  author  has  been  employed  during  the  past  year 
and  a  half  by  the  Organic  Chemicals  Division  of  The  Glidden  Company 
and  it  was  a  happy  circumstance  that  they  were  interested  In  using 
computer  control  of  a  crude  tall  oil  distillation  process  at  their 
Port  St.  Joe,  Florida  plant.   The  work  undertaken  was  primarily  of 
interest  for  its  theoretical  aspect  but  it  was  most  gratifying  to  be 
able  to  actually  test  the  theory  in  a  real  world  situation. 

Specifically,  the  problem  discussed  in  this  work  is: 

1.  To  describe  in  mathematical  terms  the  multicomponent 
distillation  process  as  it  applies  to  the  tower  designated  A-30. 

2.  To  develop  a  calculation  procedure  which  will  allow  the 
mathematical  model  to  be  solved. 

3.  To  fit  the  mathematical  model's  predicated  composition 
profiles  to  those  obtained  experimentally  from  the  distillation  column. 

4.  To  use  the  mathematical  model  to  develop  insight  into  the 
problem  of  digital  computer  control  of  the  process. 
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5.  To  develop  and  implement  the  mathematical  model  to  provide 
predictive  control  of  the  actual  process  located  at  a  distance  of  about 
300  miles  from  the  digital  computer. 


CHAPTER  II 

MATHEMATICAL  PROCEDURES 

A  schematic  diagram  of  the  column  (designated  A-30)   is  shown 
in  Figure  1.   The  operation  taking  place  in  A-30  consists  of  the  low 
pressure  distillation  of  a  mixture  of  fatty  acids  and  organic  impurities 
known  as  unsaponif iables  (unsaps).   The  feed  to  A-30  is  the  fatty  acid 
cut  of  crude  tall  oil  consisting  in  the  most  part  (807.,)  of  oleic  and 
linoleic  (C-18)  acids.   There  is  also  an  appreciable  amount  of  palmitic 
(C-16)  acid  together  with  higher  and  lower  straight  chain  acids  and 
the  unsaps.   The  commercial  object  of  the  process  is  to  produce  a  fatty 
acid  with  low  unsaps  (approx.  1%)  concentration  which  is  essentially 
free  of  palmitic  acid.   Steam  is  introduced  in  the  reboiler  to  allow 
the  distillation  to  take  place  at  a  low  temperature  by  reducing  the  par- 
tial  pressure  of  the  organics  and  thereby  saving  the  product  from 
thermal  degradation.   The  feed  is  heated  by  a  feed  preheater  before 
entering  the  column  and  the  overhead  vapors  are  condensed  in  an  internal 
type  condenser. 

A.   The  Mathematical  Model 
The  mathematical  model  is  derived  by  using  the  method  of 
Sorel  (19)  in  applying  energy  and  material  balances  around  separate 
sections  of  the  column. 

1.   Assumptions 

The  following  assumptions  are  made  in  the  basic  model: 

1.  Steam  is  always  a  superheated  vapor. 

2.  The  liquid  phase  and  the  vapor  phase  form  ideal  solutions 
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Figure  1,   Schematic  Diagraiii 
of  A-30  Distillation  Column 
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(i.e.  no  heat  of  mixing)  so  that  the  total  enthalpy  of  i  li  uid  or  vapor 
stream  may  be  treated  js  an  additive  function  o£  pure  component  enthil- 
pies  . 

3.  The  e-iuilibrium  may  be  represented  by  Raoult  s  law. 

4.  The  heat  capacity  and  heat  of  vaporization  for  each  component 
are  functions  of  temperature  only. 

5.  Tlie  efficiency  of  each  component  is  equal  and  constant. 

6.  There  is  no  heat  loss  to  the  atmosphere. 

7.  The  feed  is  a  saturated  liquid. 

8.  The  liiuid  and  vapor  streams  within  the  column  are  at 
their  bubble  and  dew  points,  respectively. 

2_; Derivation  of  the  Basic  Equations 

Taking  a  total  material  balance  around  the  reboiler  (envelope  T 
of  Figure  1)  we  obtain 

L_  +  S   =  B  +  I'  (1) 

2    o        1 

The  individual  component  balances  for  the  organics  can  be 

written  as 

^2^2  =  ^^1  +  ^^1  ^'^ 

and  for  water  vapor 

S   =  y  y   =  S,  (3) 

o    wl      1 

The  energy  balance  for  the  reboiler  is 

%   "■  '\H   "-   "so^o  =  \'   +  »1^  ^""^ 

The  total  enthalpy  of  any  liquid  stream  may  be  represented  bv 
c 

h   =  Vh.  X.  (5) 

n    ^  m  m 
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and  for  any  vapor  stream 


H  =  \  H.  y.   +  H 


n 

1=1 


(6) 


where  c  is  the  total  number  of  organic  components, 


Combining  Equations  4,  5,  and  6,  we  obtain 
c  c 

i=l  i=l         i=l 


+  >  h.„x.-L^  +H   S   =  \  h.,x.  R  +  )  H.,h.,V- 
^  12  i2  2    so  o    ^  ^1  il     [__,    il  ^1  1 


(7) 


Substituting  the  c  equations  represented  by  2  into  7  and 

rearranging  gives 

Q^  +  S  (H   -  H  ,)  -  BV(h    -  h._)x., 
S    o   so     si      ^_^       il     i2   il 

Vj  =  5 ^— (8) 


X^\l-  ^2>>^il 


1=1 

If  the  temperature  in  the  reboller  is  known,  we  may  calculate 
the  composition  of  the  vapor  leaving  the  reboller 

p. ,x. ^ 

The  pressure  in  the  reboller  may  be  calculated  by  summing 

Equation  9  over  c  and  using  Equation  3 

c 

)  p  ^1 

c       .^ 

^='        ^  (10) 


)/ll  =   ^    =  ^  ■  ^w 

i=i 

c 
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calculated    from  Equation  8 

c  c  ,v 


'o4(Hil    -    h2^ 


H    ,)    -    B V  (h.,    -    h.„)x.. 

si  ;       ii        i2ij 


(12) 


i  =  l 


Solving  for  tt,  in  Equation  12,  we  get 
c 

Y  -     -    ^^■ 


TT,  =  )  p  X  + c-:: :: (i3) 

^        (J   +  S  (H   -  H  ,)  -  BV(h.T  -  h.„)x.^ 
i=l         S    o   so    si      /    il    i2   il 

All  the  quantities  on  the  right  hand  side  of  Equation  13  are 
known  or  assumed  so  we  may  calculate  the  reboiler  pressure. 
The  order  of  solution  for  these  equations  is: 
1.    Obtain  tt,  from  Equation  13 

i' 

3.  Obtain  V   from  Equation  8 

4.  Obtain  L„  from  Equation  1 

5.  Obtain  x.„  from  Equation  2. 

Taking  energy  and  material  balances  around  the  bottom  of  the 
column,  below  the  side  stream  (envelope  II  in  Figure  1)  and  making  use 
of  the  equilibrium  and  Murphree  efficiency  relationships, 

y!   =  p'!'  X.  /tt  (9) 

^  m    in  xn     n 

y.  =  y.   .,  +  E„.  (y.  -  y.   ,)  (14) 

■^in   ^in-1    Vm  "^  m   -'m-l 

we  may  derive  exactly  the  same  sort  of  relationships  as  for  the  reboiler 

Defining  some  intermediate  quantities  as: 
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c 

V  * 

C,  =  )  E„.  p.  X.                                    (15) 

1    i_^  Vm  in  m 

1=1 


\S' 


h,  _^,)E,,,  p.  X.  (16) 

in    in+1   Vin*^in  m 


i=l 
c 

^4=  X(l-^in>>'in-l  (1«> 

1=1 
c 

C^  =  V(H.   -  h.  ^,)(1  -  E„.  )v.   T  (19) 

5    zl   1"    m+l       Vm  •  in- 1 

i=l 


We  may  write: 


(20) 


7T^  =  (C^  +  SC2/D^)/(1  -  C^  -  S.C^/op  (21) 

c 

V  =  D^/  y(H.   -  h.  _^-)y.  (22) 

n     1  [j   in    in+1   m 

i=l 
^,  =  B  +  V  -  S  (23) 


^n+1  =  (^-il+Vin)/Vl  (24) 

In  the  same  manner  as  above,  we  can  derive  the  equations  which 
apply  between  the   sidestream  and  the  feed  plate. 

Defining 

V  -       - 
C,  =  )  (H.  ^.  -  h.  )x,  (25) 

6    ^   in+l    ic   ic 

1=1 

°9  =  ^Q  +  (".^  -  H.  )S  +  BC-  +  CFA  -  C.  (26) 

2    S     so    sn       3  6 
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We  may  write: 


■n^  =   (C^   +  SC2/D2)/  (1   -    C^   -   SC5/D2)  (27) 


'n  =  ^^2  /  X    ^»in   -   ^n+l^^in 
i=l 


(28) 


L    .,    =  B   +  V      +  CFA   -    S„  (29) 

n+1  n  o 

The  above  equations  (27-30)  are  used  for  n  <f . 
Taking  energy  and  material  balances  around  the  entire  column, 
we  can  calculate  a  value  for  the  condenser  heat  duty. 


t  i;  t. 

^c  =  ^Z  ^iF^F  -  (^X^il^l  ^  ^^"  PiCFA^i 


iCFA 
i=l  i=l  i=l 

c 

+  DYh.^x._)  +  Q^  +  (H  „  -  H   )S    .  (31) 

/   iD  iD      S      SO     sD  o 

i=l 
By  taking  balances  around  the  condenser  system,  the  vapor 

composition  leaving  the  top  plate  may  be  calculated: 


^  i=l 

Vl  =  \  +  D  4-  S  (33) 

The  pressure  at  the  plate  can  be  calculated  by 


i^  P 


iD-1 
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The  following  equations  are  used  to  calculate  conditions  in  the 
top  of  the  column: 

\n   =  Vin^PL  ^^^^ 

^n  =  ^n+1  -  \in(^n+l  "  ^n>  ^^'^ 

c 


r^^  +  (H    -  H    .)S  +  D  ^   (H,_  -  H    i)x.„ 
MC     so    sn-1       j^^        ID    in-1   iDJ 


(38) 


V  (H,   ,  -  h.  )x. 
/   ^  in-1    m   m 


1=1 


V    =  L  +  D  +  S  (39) 

n-1    n 

y.   ,  =  (L  X.   +  Dx.^)/V   ,  (AO) 

■^in-1     n  m     lD   n-1 


vi-^y^  <"' 


Equations  36  and  37  are  repeated  by  incrementing  n  and  38-41  are 
repeated  until  n  =  f  +  1.   At  the  feed  plate,  an  average  composition 
of  the  liquid  entering  the  feed  plate  is  calculated  by  Equation  42, 

The  liquid  which  would  be  in  equilibrium  with  the  vapor  leaving 
the  feed  plate  is  calculated  by  Equation  36  with  n  =  f  and  the  liquid 

leaving  the  feed  plate  as  calculated  from  the  top  la 


r  -  =  X.  .   -  E,  _(x,  .   -  x*^)  (43) 

If     imix    Lif   imix    if 
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Wlien  the  compositions  calculated  from  Equation  A3  and  Equation 
30  are  equal,  the  correct  top  and  bottom  compositions  have  been  selected 
and  the  solution  has  converged. 

B.   Calculation  and  Convergence  Techniques 
One  of  the  basic  problems  in  any  multicomponent  distillation 
calculation  is  that  of  convergence  to  the  proper  solution.   This 
problem  is  especially  acute  if  a  modern  digital  computer  is  used,  since 
manual  intervention  is  for  all  practical  purposes  impossible  without 
great  losses  in  time  and  money.   In  any  hand-calculated  solution,  much 
subconscious  judgment  comes  into  play  allowing  the  calculator  to  abandon 
fruitless  paths  without  carrying  the  solution  very  far.   However,  a 
computer,  no  matter  how  fast  or  sophisticated,  cannot  exercise  the 
engineering  judgment  of  the  most  unskilled  practitioner  and  therefore 
must  be  programmed  in  a  foolproof  method  which  assures  convergence  to 
a  proper  solution.   By  convergence,  it  is  meant  that  by  taking  sufficient 
trials,  the  error  may  be  reduced  to  within  any  predetermined  limits, 

1.   Previous  Work 

Smith  (9)  has  an  excellent  summary  of  the  various  methods  of 
choosing  new  variables  for  the  distillation  calculations  and  two  books 
(10,24)  have  recently  been  ^^/ritten  on  this  subject  alone. 

Perhaps  the  oldest  method  in  use  is  that  of  direct-iteration 
which  uses  the  results  of  the  previous  trial  without  any  manipulation 
to  start  the  next  trial.   Lewis  and  Matheson  (20)  in  their  original 
paper  illustrate  the  direct  iteration  technique  as  do  Thiele  and 
Geddes  (23). 

Bonner  (22)  was  one  of  the  first  to  present  a  computer  program 
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for  the  solution  of  multicomponent  distillation  problems.   He  programmed 
the  basic  Lewis-Matheson  prooedure  with  slight  adjustment  for  machine 
computation.   However,  his  technique  for  convergence  was  direct- 
iteration. 

Lyster,  et  al.  (33)  and  Holland  (24)  have  presented  the  Theta 
method  of  convergence  and  used  it  with  the  Thiele-Geddes  method  of 
calculation.   Since  many  articles  have  appeared  on  this  technique,  a 
brief  description  of  it  seems  appropriate. 

Let  us  relate  the  assumed  (trial)  and  correct  values  of  product 
ratios  as  follows:   For  any  component, 


a^   -^(^) 


(44) 
correct         "   calculated 

From  the  individual  component  material  balances  and  the  require- 
ment that  the  overhead  flow  rate  is  fixed,  we  may  write 


D  =  y  (d.) 

^    1  correc 


(45) 


(46) 


i=l 

If  the  final  and  assumed  values  were  known,  we  could  calculate 
the  set  of  6.'s   which  satisfy  Equation  44;  however,  this  is  not  the 
case.   We  can  find  a  single  6   which  will  satisfy  Equations  45  and  46, 
if  the  subscript  i  is  dropped  from  the  9   in  Equation  44.   This  single 
9   can  be  considered  an  approximation  of  each  value  in  the  set  of  9.'s. 
It  is  argued  that  this  is  a  reasonable  approach  because  all  the  9.'s  must 
equal  one  when  the  assumed  value  is  the  correct  value. 
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Solving  E4uations  44-46  and  eliminating  the  correct  values  of 

products  (the  unknowns)  ^we  can  obtain  E-juation  47. 

c           p 
g(9)  =  V  TTT^ D  =  0  (47) 


V  i  ■'calculat 


i'=i  ^        .  ^    . 

\     i    ■'calculated 

The  desired  value  of  0   is  the  positive  root  which  makes  g(0)  =  0, 
Once  this  6   is  found  (Holland  gives  several  standard  methods  for  finding 
the  roots  of  algebraic  equations)  the  corrected  values  of  d.  and  b  may 
be  calculated  from  Equations  48  and  49. 

(d.)      ^  ,= 7-5T-T (48) 

1  corrected    1+0 

(b.)        .  =  0  ^T^  ^  (d.)      ,  ,        (49) 

^  1  corrected       d.    J     ^      i  .  j   I'corrected 
1  /calculated 

The  Thiele-Geddes  iteration  is  then  repeated  until  convergence 
is  obtained.   For  a  more  complete  discussion  of  this  method  and  its 
usefulness,  the  reader  is  referred  to  the  original  papers  or  Holland's 
book  (24)  which  deals  with  the  Theta  method  in  great  detail. 

Greenstadt,  Bard  and  Morse  (34)  developed  a  very  complex  IBM  704 
computer  program  to  solve  the  multicomponent  distillation  problem. 
Conditions  at  the  top  and  bottom  of  the  column  are  assumed.   Calcula- 
tions are  carried  out  until  the  "feed"  or  "match"  plate  is  reached.   If 
compositions  are  not  the  same,  as  calculated  from  the  top  and  bottom, 
an  "appropriately  framed  form  of  Newton's  method"  is  used  to  make  this 
connection.   The  calculations  are  then  repeated  until  matching  conditions 
are  satisfied.   This  program,  while  very  sophisticated  from  a  machine 
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point  of  view,  is  rather  obscure  for  the  average  chemical  engineer.   It 
uses  a  machine  type  language  rather  than  the  more  readily  understood 
programming  languages,  such  as  FORTRAN,*  and  hence  the  chemical  engineer- 
ing clarity  is  sacrificed  for  the  sake  of  machine  considerations. 

Deland  and  Wolf  (35)  present  a  procedure  for  solving  distilla- 
tion problems  on  an  electronic  analog  computer  by  replacing  the  normal 
algebraic  equations  with  differential  equations  and  obtaining  the  steady- 
state  solution.   A  direct-iteration  procedure  is  used  to  correct  the 
calculated  compositions  by  varying  feed  rates  and  number  of  plates  in 
the  column. 

Peiser  (36)  presents  a  convergence  technique  which  he  claims 
is  rapid,  with  only  the  most  difficult  problems  taking  seven  or  eight 
iterations.   He  considers  the  distillation  column  to  be  made  up  of  two 
lumped  units,  rectifying  and  stripping  sections.   Consider  Figure  2 


Fx, 


Figure  2.   Block  Diagram 
of  a  Distillation  Tower 
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The  individual  componenL  material  balances  which  must  be  satisfied 
are  shown  in  Equations  50  and  51. 

ViT  "  '^iD  =  ViB  -^  ^-^G  <50) 

^^T  ^  ^'^M  =  ^'B^'iB  -^  ^^B  <51) 

If  these  equations  are  not  satisfied,  then  new  trial  values  of  B,  x.„. 

Id 

D  and  x   must  be  picked.   This  is  done  by  treating  the  feed  as  if  it 
was  flashed  and  we  seek  a  relationship  of  the  form 


SX 


(52) 


where  the  7.'s  are  so  defined  as  to  equal  one  when  Equations  50  and  51 
are  satisfied.   Alpha  is  the  relative  volatility  of  component  i.   After 
a  certain  amount  of  admitted  heuristic  argument,  Peiser  defines  the 


rViB_l^        VlB    -^  ^^B  1 
i      LVir^^^D       Vir+\iMJ 


T^^^ile  this  form  is  quite  arbitrary,  Peiser  claims  rapid  con- 
vergence in  all  cases. 

Hansen  (10)  presents  two  methods  of  convergence  in  his  book 
of  computer  programs.   His  Method  I  is  the  same  as  Holland's  (24) 
Theta  method  and  his  Method  II  is  derived  from  the  basic  differential 
equations  of  the  dynamic  system. 

Hansen  uses  the  following  equation 

, ,    V   ,  K   ,x   ,  +  L  .  ,x  ,, 
/   \i^+l    n-1  n-1  n-1    n+1  n+1  ,,,. 

(   )    . (54) 
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to  calculate  the  value  of  x   for  the  r+lth  trial.   An  additional 
n 

assumption  must  be  made  to  start  this  solution  and  that  is  the 
composition  at  every  stage.   Hansen  has  found  it  satisfactory  to  use 
feed  liquid  at  its  bubble  point  on  each  plate  as  a  first  assumption. 
Equation  54  has  the  exact  form  of  the  well-known  flash  equation,  so 
Hansen  has  termed  this  method  as  the  "Method  of  Successive  Flashes." 

Other  convergence  techniques  include  those  of  Bonner  (22) 
and  Ball  (28). 

2.   The  Power  Ratio  Convergence  Technique 

A  method  developed  in  this  work  which  is  presented  here  makes 
use  of  a  ratio  of  compositions  at  the  feed  plate  as  calculated  from  the 
top  and  bottom  by  a  plate-to-plate  calculation.   This  ratio  raised  to  a 
power  which  acts  as  a  "gain  factor"  in  eliminating  the  error  between 
top  and  bottom  calculations  has  been  found  to  be  particularly  useful  in 
this  work  and  appears  to  have  a  considerable  potential  for  wide-spread 
application  in  multicomponent  distillation  work.   This  method  was 
developed  because  of  difficulties  which  were  encountered  in  trying  to 
use  Holland's  Theta  method  on  A-30  tox>;er .   The  difficulties  hinged  on 
the  fact  that  A-30  is  a  complex  tower  having  restrictions  on  the  sidestream 
composition  and  having  an  inert  component  (steam)  passing  through  the 
whole  tower.   Insufficient  work  was  done  to  prove  that  Holland's  method 
fails  or  must  be  modified  for  this  tower,  but  it  was  clearly  very 
difficult  to  implement.   It  is  also  clear  from  Holland's  writings  (24) 
that  complex  columns  require  special  arrangements  of  his  equations  that 
often  depend  on  the  particular  column  being  considered. 
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Only  one  reference  has  been  found  (37)  in  which  the  use  of  a 
ratio  a6  the  basis  for  a  convergence  technique  is  mentioned.   The  ratio 
is  used  without  the  gain  factor  and  in  an  entirely  different  fashion 
from  that  proposed  here.   There  is  no  report  of  its  use  except  with  a 
single  simple  tower. 

In  the  Power  Ratio  Convergence  Technique,  the  form  of  the 
distillation  equations  is  unimportant;  all  that  is  necessary  is  that 
one  be  able  to  calculate  compositions  from  both  directions,  starting 
with  composition  values  in  the  top  and  bottom  streams.   The  choice  of 
the  particular  stream  on  which  the  convergence  is  made  is  unrestricted, 
as  the  method  works  equally  well  on  vapor  streams  or  other  lijUid  streams, 
The  feed  plate  is  a  convenient  choice  because  this  avoids  the  possibility 
of  computing  negative  values  from  the  material  balance  equations.   Wien 
we  compute  negative  flow  rates  or  composition,  we  move  mathematically 
into  a  physically  unattainable  region  which  violates  the  implied 
assumption  in  all  distillation  problems  of  a  real-world  situation  with 
positive  flow  rates  and  compositions. 

If  we  designate  the  liquid  leaving  the  feed  plate  as  x^^.^,  when 

calculated  from  the  top  and  as  x.^„  when  calculated  from  the  bottom  we 

itB 

can  compute  a  ratio  for  each  component  i, 

y.    =  ^^  (55) 

'        ^fB 

This  ratio  is  equal  to  1  when  both  computed  compositions  are  the  same. 
If  we  examine  the  ratio  relative  to  calculated  composition  profiles,  we 
find  that  it  is  greater  than  1  when  the  profiles  do  not  meet,  as  shown 
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in  Figure  3,  and  less  than  1  when  the  profiles  overlap  one  another,  as 

shown  in  Figure  4.   A  ratio  greater  than  1  indicates  that  the  guess  of 

x.„  (bottom  composition)  was  too  small  and  x.„  (top  composition)  vvas 
lis  il 

too  large,  while  a  ratio  less  than  1  indicates  the  exact  opposite. 

The  correction  to  make  these  profiles  meet  at  the  feed  plate  can 

be  shown  if  we  consider  the  bottom  half  of  the  column  alone.   When  the 

column  profile  is  plotted  with  composition  on  a  logarithmic  scale  and 

stage  number  on  an  arithmetic  scale,  for  different  trial  values  of  x . ^ , 

IB 

an  important  property  of  the  composition-position  function  is  observed. 
The  second  profile  appears  to  be  a  linear  displacement  of  the  first  with 
respect  to  the  log  axis.   As  shown  in  Figure  5,  a  change  in  x.„  produces 

Id 

an  approximately  equal  change  in  x.^„. 

irrJ 

We  can  express  the  column  profile  for  any  trial  as  a  function 
of  log  X. 

N  =  0(log  X.)  (56) 

and  transform  this  function  into  the  profile  for  any  other  trial  by  the 
linear  translation 

log  x^  =  log  x_.  +  log  a  =  log  ax.  (57) 

The  factor  a  is  a  function  of  the  ratio  term  defined  by 
Equation  55  since  the  ratio  is  the  percentage  on  the  log  axis  by  which 
we  want  to  displace  the  profile.   The  ratio  is  the  factor  by  which  all 
points  in  the  bottom  profile  must  be  multiplied  in  order  to  displace  the 
feed  plate  composition  to  meet  the  previous  value  of  the  feed  plate 
composition  which  was  calculated  for  the  top  profile. 

The  top  and  bottom  compositions  are  directly  related  by  the 
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overall  material  balance  so  any  change  in  one  effects  the  other. 

Therefore,  if  we  shift  the  bottom  profile  to  meet  the  top  profile,  the 

top  profile  also  moves  and  the  profiles  could  possibly  pass  one  another 

on  the  next  trial,  especially  if  the  top  and  bottom  compositions  are 

of  the  same  order  of  magnitude.   To  allow  for  this  problem,  each  ratio 

is  raised  to  a  power,  the  value  of  which  is  selected  based  on  past 

history  of  the  ratio,  that  is  on  whether  the  ratio  is  oscillatory  about 

1  or  not.   We,  therefore,  define  a  gain  factor  g.  by  the  equation 

^i 
k.  =7^  (58) 

in  order  to  control  the  corrective  action. 

Since  x.^,  and  x,^  are  related  by  tl 
Id       1  i 

some  method  must  be  devised  to  allow  both  the  restrictions  of  match  and 
component  material  balance  to  be  satisfied.   One  obvious  method  is  to 
apply  all  the  ratio  corrections  to  either  the  top  or  bottom  stream  and 
calculate  the  other  stream  composition  using  the  overall  component 
material  balance.   Using  this  technique  has  the  disadvantage  that  a 
component  whose  concentration  is  high  in  the  stream  corrected  by  the 
ratio  may  be  changed  so  that  the  material  balance  gives  a  negative 
concentration  in  the  other  stream.   This  difficulty  is  avoided  by 
correcting  those  components  which  concentrate  in  the  top  (lights)  in 
the  bottom  and  those  components  which  concentrate  in  the  bottom  (heavies) 
in  the  top.   In  addition,  some  consideration  should  be  made  of  the  side- 
stream  concentrations.   For  a  sidestream  which  is  located  in  the  bottom 
section  of  the  tower,  it  appears  logical  to  apply  the  same  correction 
as  we  have  made  to  the  bottom  stream,  while  if  the  sidestream  is 
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located  in  the  top  section  of  the  tower,  we  should  make  the  same 
correction  as  we  have  made  in  the  top  stream. 

Based  on  these  considerations,  the  following  equations  were 
used  to  correct  the  assumed  top  and  bottom  compositions  for  A- 30 
column  which  had  one  sidestream  located  in  the  bottom  section  of  the 
column . 
For  the  light  components 

^B  =  ^B^  <59) 


and  by  material  balance, 


N    ^^F  ■  ^^^^'^JcFa/'^Jb  -^  «>^B 


"iD 
For  the  heavy  components 


(60) 


iT    iT   i 

and  by  material  balance, 

M       Fx._  -  Dx.^ 
N         3.F iD 


^  -x?/.  (61) 


^^^^^cfa/^b  ^  ^^ 


(62) 


In  the  event  that  Equations  60  or  62  give  a  negative  number, 
the  previous  value  of  that  variable  is  used  in  the  next  trial. 

It  has  been  found  that  the  rate  of  convergence  of  this  technique 
is  a  function  of  the  gain  factor.   At  a  fixed  gain  of  1  for  all  compo- 
nents (applying  the  entire  ratio  correction)  the  method  is  oscillatory; 
therefore,  to  begin  the  solution,  gains  of  less  than  one  are  used  with 
the  ratio  terms  being  monitored  in  the  computer  program.   Any  repeated 
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change  in  ratio  from  greater  than  1  to  less  than  1  (which  indicates 
oscillation)  causes  a  reduction  in  the  size  of  the  gain  factor. 
However,  it  is  desirable  to  keep  the  gain  as  high  as  possible  in  order 
to  reduce  the  number  of  column  calculations  necessary  for  convergence; 
therefore,  the  gain  is  increased  on  those  components  which  are  not 
oscillatory. 

Figures  6  through  8  show  the  effect  of  this  gain  factor  upon 
convergence  time.   These  plots  are  of  the  ratio  for  same  three 
components  vs.  iteration  number.   Figure  6  shows  a  gain  of  0.5,  while 
Figure  7  illustrates  the  oscillatory  behavior  with  a  constant  gain  on 
each  component  of  unity.   The  data  which  are  plotted  in  Figure  8  was 
made  with  a  constant  gain  equal  to  0.9.   It  is  interesting  to  note  that 
the  reduction  in  error  is  almost  twice  as  much  per  iteration  in  this 
case  as  it  is  with  a  gain  of  0.5. 
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Figure  6.   Ratio  with  Constant  Gain  of  0.5 
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Figure  7.   Ratio  with  Constant  Gain  of  1.0 
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CHAPTER  III 
MODEL  FITTING  AND  EXPERIMENTAL  RESULTS 
To  be  of  practical  use,  a  theoretical  model  must  describe  the 
process  with  reasonable  accuracy.   One  of  the  identifying  attributes  of 
most  industrial  operations  is  the  lack  of  accurate  and  reliable 
physical  data.   This  lack  of  data  must  be  overcome,  at  least  to  the 
extent  of  having  the  model  duplicate  the  actual  plant  process  -within 
acceptable  limits.   This  was  done  in  this  work  by  varying  the  constants 
in  the  vapor  pressure  relationships  until  a  good  fit  of  the  experimental 
plant  data  was  obtained.   These  constants  were  then  verified  by  an 
independent  set  of  plant  data. 

A .   Physical  Data 


An  examination  of  the  model  as  derived  in  Chapter  II  will  show 
that  there  are  several  parameters  which  can  be  changed  to  obtain  a 
good  fit  for  the  model.   The  two  most  important  parameters  used  in  the 
model  are  the  vapor  pressure  of  each  component  and  the  plate  efficiencies, 
Therefore,  it  became  necessary  to  define  the  components  in  the  tall  oil 
fatty  acid  and  determine  or  assign  vapor  pressures  to  each  of  them. 

The  feed  to  A-30  column  consists  of  the  fatty  acid  (FA)  portion 
of  crude  tall  oil  (CTO) .   The  FA  is  composed  of  the  straight  chain  fatty 
acids,  both  saturated  and  unsaturated,  with  carbon  numbers  of  14  and 
higher  and  certain  structurally  unknown  substances  designated  as  unsaps . 
From  vapor  phase  chromatography  (VPC)  the  feed  stock  to  A-30  was  split 
into  4  acids  and  5  unsaps.   The  acids  are  designated  as  C- 14  (Myristic), 
C-16  (Palmitic),  C-18  (Oleic),  and  C-20  (Sterric)  and  the  unsaps  as 
L-1  through  L-5  with  molecular  weights  assigned  as  239,  290,  312, 
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318,  and  332,  respectively.   The  VPC  method,  of  determining  the  compo- 
sition of  FA  presented  in  Appendix  D,  was  developed  in  the  Port  St.  Joe 
Laboratory  and  is  faster  and  more  detailed  than  the  standard  method  (38) 

Vapor  pressure  data  on  pure  components  were  collected  from  as 
many  sources  as  possible  and  are  shown  in  Appendix  B.   These  data 

were  fit  to  an  Antoine  type  equation: 

A. 
log  p."  =  7^-^  +  B.  (63) 

'i 

By  the  use  of  the  Clausius-Clapeyron  equations,  a  relationship 

was  derived  to  give  the  latent  heat  of  vaporization  on  a  function  of 

temperature  (39). 


1      log  c  10   1  I   t  +  C.  J 


(64) 


Liquid  heat  capacity  data  were  obtained  from  the  literature 
for  those  substances  which  were  known  (straight  chain  acids)  and  an 
average  c  .  based  on  the  fatty  acids  were  used  for  the  unsaps .   All 
of  the  data  obtained  on  these  substances  could  be  reasonably  related  to 
temperatures  by  a  linear  relationship. 

c  .  =  a.  +  b.t  (65) 

pi      L      1 

Saturated  liquid  enthalpies  were  calculated  by  Equation  66  and 
vapor  enthalpies  were  calculated  by  Equation  £7 . 
t 

-      r " 

h.   =   /   (a.  +  b.t)dt  (66) 

H.   =h,   +X.,  (67) 

m    in     itn 

Using  the  available  literature  values  for  vapor  pressure  and 
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assuming  100  per  cent  efficiency,  a  steady-state  calculation  on  the  com- 
puter was  made  with  flow  rates,  feed  composition,  and  reboiler  heat  duty 
corresponding  to  an  actual  experimental  run  on  the  column.   The 
efficiencies  were  lowered  until  the  trends  of  the  computed  profiles 
tended  to  be  the  same  as  the  observed  profiles.   The  constants  in  the 
vapor  pressure  curves  were  then  varied  until  a  reasonable  fit  of  the 
experimental  data  was  made.   The  final  physical  constants  used  in  this 
work  are  also  shown  in  Appendix  B. 

B.   Experimental  Work 

Two  complete  sets  of  data  were  taken  on  the  tower.   A  complete 
set  of  data  consists  of  samples  of  the  feed,  top,  side,  and  bottom 
streams   of  the  column,  8  samples  from  trays  within  the  tower,  of  flow 
rates  of  all  streams,  temperatures  at  selected  plates  in  the  column,  and 
the  reboiler  heat  duty.   One  set  of  data  (040165)  was  used  to  obtain  the 
proper  physical  constants  for  the  vapor  pressures  of  each  component  as 
described  in  Part  A  of  this  chapter.   The  other  set  of  data  (050665)  was 
used  to  verify  the  model's  fit  to  the  actual  process.   In  addition,  a  par- 
tial set  of  data  (i.e.,  everything  but  internal  plate  samples)  was  obtained 
at  rates  which  differed  drastically  from  the  normal  operating  conditions. 

Table  1  shows  the  experimental  analysis  of  the  samples  collected 
for  Run  040165.   Table  2  shows  the  operating  conditions  for  Run  040165. 
Figure  9  shows  the  total  unsaps  concentration  as  calculated  by  the 
model  and  as  obtained  experimentally  for  Run  040165.   Figures  10-18 
show  each  individual  component  profile  for  Run  040165  as  obtained  experi- 
mentally and  as  calculated  with  the  mathematical  model. 

The  disagreement  between  the  experimental  points  and  the 
calculated  curves  for  L-1  (Figure  11)  and  L-2  (Figure  13)  is  probably 
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TABLE  2 
Operating  Condi tions^-' 
Run  040165 

Feed  rate  100  units/hr, 


Bottoms 


Steam 


5.2  units/hr. 


Sidestream  70.8  units/hr, 

Heads  **  21.8  units/hr, 


39 .4  units/hr. 


Reboiler  Heat  Duty  8.67  x  10   Btu/hr. 


Reflux  Temperature 


85^C 


*  Flow  rates  are  based  on  100  units/hr.  of  feed 

**  Overhead  product  is  designated  as  Heads  throughout 
this  work. 
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Figure  14.   C-18  Profile  Run  040165 
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due  to  the  analysis  technique.   The  feed  to  the  tower  actually 
contains  a  very  small  amount  rosin  acid   (1-27,)  which  is  left  after  the 
split  between  rosin  acid  and  FA  is  made  on  the  CTO  in  a  tower  immediately 
upstream  of  A-30  tower.   A  rosin  acid  determination  on  the  bottoms  of 
the  column  shows  that  almost  all  of  these  acids  go  out  of  the  bottom  of 
the  tower  and  that  none  can  be  determined  above  plate  8  (AO) .   In 
addition,  the  retention  position  occupied  by  these  unsaps  on  a  VPC  is 
the  same  as  that  occupied  by  the  rosin  acids  present  in  CTO.   Therefore, 
the  reported  unsaps  concentrations  in  the  bottom  of  the  tower  for  the 
light  unsaps  is  most  likely  in  error  and  what  really  are  present  are 
rosin  acids.   The  disagreement  shown  in  the  C- 14  profile  (Figure  10) 
in  the  bottom  of  the  column  probably  comes  from   polymerizations  of  the 
C-14's  in  the  bottom  due  to  high  temperature  and  long  holding  times  and 
a  cracking  which  occurs  during  the  VPC  analysis. 

The  agreement  between  experimental  points  and  calculated 
curves  for  the  other  components  is  within  the  experimental  error  of 
the  analysis  technique  and  is  rather  good  for  actual  plant  process  data. 

Table  3  gives  the  experimental  compositions  data  for 
Run  050665  and  Table  4  gives  the  operating  conditions  for  this  same 
run.   Figure  19  shows  the  total  uns.ip  profile  and  Figures  20-28  show 
the  individual  component  profiles  for  Run  050665. 

The  agreement  between  calculated  and  experimental  points  is 
better  for  this  run  for  L- 1  (Figure  21)  and  L-2  (Figure  23)  than  for 
the  previous  run.   Except  for  the  bottom  stream,  no  L- 1  is  detected 
below  tray  8  and  the  probable  reason  is  that  the  feed  stock  in  this 
case  was  exceptionally  low  in  rosin  acids  (40). 
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TABLE  4 
Operating  Conditions* 
Run  050665 


Feed  rate 

Bottoms 

Sidestream 

Heads 

Steam 

Reboiler  Heat  Duty 

Reflux  Temperature 


100  units/hr. 
4.7  units/hr. 
75.0  units/hr. 
20.2  units/hr. 
39,6  units/hr. 
7.16  X  10^  Btu/hr, 


*  Flow  rates  are  based  on  100  units/hr.  of  feed 
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Figure  19.   Total  Unsaps  Profile  Run  050665 
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Figure  24.   C-18  Profile  Run  050665 
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Run  060365  was  made  at  one-half  the  rates  normally  used  on 
the  column.   Table  5  shows  the  operating  conditions  for  this  run.   The 
mathematical  model  predicted  an  unsaps  concentration  in  the  sidestream 
of  0.72  mole  per  cent  while  the  actual  experimental  value  of  unsaps  was 
0.67  mole  per  cent.   This  was  a  rather  drastic  test  for  the  model, 
since  the  flow  rates  for  this  run  were  50  per  cent  of  those  in  the  run 
used  to  fit  the  model.   In  addition,  the  control  model,  as  will  be 
discussed  in  Chapter  IV,  predicted  a  heads  rate  of  17.8  units/hour 
to  bring  the  unsaps  concentration  to  1.4  per  cent.   At  the  time  of  this 
run,  the  computer  and  the  associated  equipment  necessary  for  routine 
transmission  of  data  between  Port  St.  Joe  and  Jacksonville  were  not 
installed  and,  therefore,  this  control  model  was  run  after  the  fact. 
However,  the  operator,  on  the  basis  of  laboratory  reports  of  unsaps 
concentrations,  brought  the  sidestream  on  specifications  at  1.4  per  cent 
by  using  a  heads  rate  of  18.5  units/hour. 

C.   Evaluation  of  Operating  Variable  Accuracy 
One  interesting  result  from  this  work  was  an  investigation  of 
the  effect  of  accurate  values  for  the  various  operating  variables  used 
in  this  particular  model.   Most  of  the  pertinent  parameters  were  varied 
to  determine  their  individual  effects.   Since  one  of  the  important 
computed  variables  is  total  unsaps  in  the  sidestream,  the  following 
discussion  is  on  the  effect  of  the  varied  parameters  on  this  concentra- 
tion. 

The  concentrations  of  unsaps  in  the  sidestream  was  1.40  mole 
per  cent  at  an  arbitrary  set  of  standard  operating  conditions.   A  change 
in  reboiler  heat  duty  of  ±10  per  cent  gave  a  change  in  this  concentration 
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TABLE  5 


Operating  Conditions* 
Run  060365 


Feed  rate 

Bottoms 

Sidestream 

Heads 

Steam 

Reboiler  Heat  Duty 

Reflux  temperature 


100  units/hr. 
6.5  units/hr. 
60.5  units/hr. 
33.0  units/hr  . 
43.0  units/hr . 
8.05  X  10^  Btu/hr 
87°C 


*   Flow  rates  based  on  100  units/hr.  of  feed, 
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of  ±0.08  mole  per  cent  while  a  change  in  steam  flow  rate  of  ±  10  per  cent 
made  a  change  of  ±0.02  mole  per  cent.   A  change  in  feed  plate  from 
plate  16  to  plate  20  increased  the  unsaps  concentration  by  0.08  mole 
per  cent,   Kone  of  these  parameters  affected  the  distribution  of  the 
individual  components  to  any  large  degree.   The  individual  component 
profiles  were  well  within  the  experimental  error  of  the  composition 
analysis . 

Although  it  was  not  done  for  the  control  project  described 
here,  such  simulations  as  those  described  above  should  be  done  before 
completely  specifying  the  hardware  for  a  computer  control  system.   The 
results  of  these  types  of  simulation  would  indicate  which  process 
variables  are  critical  and  should  be  measured  with  a  high  degree  of 
accuracy,  which  variables  should  be  monitored  only  on  a  routine  or 
demand  basis,  and  which  variables  are  unimportant  to  the  overall  process 
and  do  not  need  to  be  monitored  by  the  computer  system. 


CHAPTER  IV 
THE  CONTROL  PROBLEM 

A.   Mathematical  Procedures 

The  mathematical  model  developed  for  A-30  Column  was  modified 
so  that  it  could  be  used  in  a  feedforward  fashion  to  predict  the 
overheads  rate  necessary  to  keep  the  unsaps  concentrations  in  the 
sidestream  within  desired  specifications.   This  modification  consisted 
of  replacing  one  degree  of  freedom  (heads  rate)  with  another  degree 
of  freedom  (unsaps  concentration)  and  obtaining  the  steady-state 
solution. 

Other  authors  (10,24)  have  suggested  that  this  problem  be 
handled  by  repeated  solutions  of  design  distillation  models  similar 
to  the  simulation  model  presented  in  Chapter  II.   Since  this  technique 
would  involve  a  multiplicity  of  steady-state  solutions  before  obtain- 
ing the  desired  operating  information,  it  seems  that  less  time  would 
be  required  if  the  calculations  converged  to  the  operating  point 
and  the  final  steady  state  simultaneously.   Therefore,  the  method 
presented  below  was  developed  to  provide  this  type  solution. 

The  ratio  technique,  in  addition  to  calculating  new  top  and 
bottom  compositions  for  the  next  trial,  also  provides  estimates  for 
the  concentration  of  any  exit  streams.   The  estimated  unsaps  concentra- 
tions in  the  sidestream  is  calculated  and  checked  to  see  if  it  is 
within  the  desired  limits.   If  it  isn't,  the  heads  rate  is  changed 

proportional  to  the  estimated  concentrations. 

unsaps 

d'^  =  D  7T1 7—  (68  ) 

specification 
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When  this  technique  is  applied  every  iteration,  the  solution 
is  very  sensitive  to  these  changes  and  the  convergence  rate  is  slowed 
down.   Figures  29  and  30  show  the  change  in  unsaps  concentrations  and 
the  heads  rate  change  for  a  calculation  using  this  technique. 

To  damp  out  the  oscillation,  the  heads  rate  correction  was 
applied  every  other  iteration  and  the  top  and  bottom  compositions  were 
recalculated  using  the  new  flow  rates  so  that  the  material  balance 
restrictions  were  satisfied.   Figures  31  and  32  show  the  unsaps 
concentrations  and  heads  rate  as  calculated  in  this  fashion. 

The  computer  listing  for  the  control  model  convergence  sub- 
routine is  shown  in  Appendix  C. 

B.   Experimental  Verification 
Once  the  simulation  of  A-30  was  made  and  the  model  described 
the  actual  column  within  experimental  limits,  the  control  model  as 
described  above  was  programmed  for  the  IBM  1620/1070  System.   There 
are  several  limitations  on  the  programming  which  are  a  result  of  the 
choice  of  this  particular  computer  system.   In  effect,  the  program 
had  to  be  "stream- lined"  in  such  a  fashion  as  to  make  it  as  fast  and 
as  efficient  is  possible.   The  first  effort  to  reduce  the  computation 
time  of  the  model  was  to  replace  the  simulation  of  "real  trays"  having 
Murphree  efficiencies  less  than  100  per  cent  by  an  equivalent  number  of 
"ideal  trays"  at  100  per  cent  efficiency.   The  main  problem  here  was  to 
maintain  a  1  to  1  correspondence  between  the  various  external  streams.   The 
model  was  reduced  to  15  ideal  plates  (vs.  29  in  the  actual  column)  with  the 
sidestream  at  plate  2  (vs.  3  in  the  actual  column).   Several  tests  were 
run  on  this  model  with  the  unsaps  concentration  and  external  flow  rates 
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identical  to  the  full  scale  model.   The  second  simplification  of  the 
model  was  to  assume  a  constant  temperature  profile  for  this  reduced 
model.   This  assumption  removed  the  necessity  of  monitoring  the  actual 
tray  temperatures  from  the  real  column  averaging  them  in  some  fashion 
for  use  on  the  reduced  column  and  then  computing  the  vapor  pressures, 
and  vapor  and  liquid  enthalpies  for  each  component  at  each  plate. 
Computer  calculations  indicated  that  slight  differences  in  temperatures, 
as  would  occur  in  normal  day  to  day  operation,  did  not  effect  the 
concentration  profile  of  the  model.   These  simplifications  of  the 
model  reduced  the  computer  running  time  by  about  40  per  cent. 

The  control  of  A-30  tower  is  initiated  in  the  following 
fashion.   A  sample  of  the  feed  stream  and  the  sidestream  is  collected 
every  4  hours  by  operating  personnel.   These  samples  are  then  taken 
to  the  laboratory  where  they  are  prepared  and  run  on  the  VPC  (see 
Appendix  D) .   The  results  from  the  VPC  are  entered  into  the  IBM  1070 
at  Port  St.  Joe.   The  computer  accepts  these  results,  calculates  the 
weight  per  cent  of  each  component  and  sends  this  information  back  to 
the  Port  St.  Joe  Plant.   The  computer  then  collects  data  from  the 
column  at  Port  St.  Joe  and  makes  a  series  of  mechanical  checks  on  key 
variables.   These  variables  are  top  column  pressure  and  temperature, 
base  column  pressure  and  temperature,  and  steam  flow  rate.   If  these 
variables  are  within  predetermined  limits,  the  first  program  for  the 
control  calculations  is  then  loaded  into  the  core  of  the  computer 
(see  Figure  33  for  the  flow  chart  of  the  A-30  control  programs).   This 
first  program  converts  the  weight  per  cent  from  the  preceding  program 
to  mole  per  cent,  compares  the  feed  analysis  and  rate  to  the  last 
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Figure  33.   A-30  Control  Programs 
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available  feed  analysis  and  rate  and  makes  the  decision  to  run  the 
control  model.   If  control  is  to  be  exercised,  the  results  are  reported 
back  to  Port  St.  Joe  as  heads  rate  calculated  in  per  cent  of  chart 
reading  for  the  use  of  the  operators.   The  heads  to  feed  ratio  is 
also  calculated  and  this  ratio  is  maintained  by  the  operators  until 
the  next  time  the  control  model  is  run.   This  is  necessary  because 
of  changes  made  in  the  A-30  feed  rate  to  control  another  column  which 
produces  the  feed  to  A-30  Column. 

Data  for  evaluation  of  the  model  and  the  control  scheme  were 
obtained  during  a  test  period  from  December  15  to  19,  1965.   In  order 
to  provide  a  basis  for  comparison,  data  vere  taken  under  operator 
control  conditions  for  the  four  days  preceding  the  actual  test.   Table  6 
shows  the  yields  for  both  of  these  periods.   These  yields  are  based 
on  flow  rates  which  are  monitored  every  40  seconds  by  the  IBM  1620/1070 
and  the  average  rate  is  computed  in  engineering  units  (pounds  per  hour) 
approximately  every  2  hours.   During  the  test  period,  there  was  a 
1  per  cent  increase  in  crude  fatty  acid  production  and  a  1  per  cent 
decrease  in  heads  production.   Based  on  current  price  differential 
between  the  value  of  crude  fatty  acid  and  heads,  this  represents  a 
savings  of  approximately  $10,000  per  month  for  the  average  crude  tall 
oil  plant  where  throughput  is  10,000  tons  per  month  (i.e.  savings  of 
$1  per  ton  of  CTO) . 

The  fatty  acid  produced  during  this  test  period  was  below  the 
maximum  unsaps  allowable  for  sale  of  this  material  and  a  further 
increase  in  yield  could  be  obtained  by  adjusting  the  unsaps  specifica- 
tions in  the  model  closer  to  the  upper  limit.   This  control  point  in 


TABLE  6 


Operating 

Yields  for  Test  Peric 

|ds 

Dates 

Units  of  A- 
Feed^' 

■30 

Units  of  CFA* 
24.24 

Units  of  Heads* 

Dec. 
11-15 

36.69 

8.59 

15-19 

36.51 

25.23 

7.51 

ised  on  total  CTO  throughput  to  plant  as  100  units 


the  model  was  kepU  low  purposely  daring  the  test  so  as  to  minimize 
the  risk  of  producing  off-grade  material. 

Table  7  gives  the  unsaps  concentrations  as  determined  by  the 
VPC  during  the  duration  of  the  test.   The  average  unsaps  concentrations 
was  1.43  per  cent  with  a  sample  variance  of  0.04.   The  control  point 
in  the  model  was  1.40  per  cent  and  a  "t"  test  showed  that  the  average 
was  not  significantly  different  from  the  control  point. 

The  control  model  was  put  on  line  successfully  at  the  Port  St.  Joe 
plant  on  a  continuing  basis  after  a  plant  shutdown  which  occurred  at 
the  end  of  December,  1965.   In  general  the  results  of  the  test  period 
have  been  confirmed  during  continuous  operation.   Improved  economic 
operation  has  been  maintained  and  operator  confidence  has  continually 
increased.   An  interesting  by-product  of  the  industrial  implementation 
oi  this  work  became  apparent  when  due  to  equipment  difficulties  (i.e.  the 
VPC)  the  control  of  the  column  was  returned  to  the  operators.   The 
method  of  control  of  the  column  (without  the  computer)  is  to  make 
adjustments  on  the  heads  rate  on  the  basis  of  one  unsaps  determination 
made  on  a  24  hour  composite  of  material.   The  operators  began  to  make 
much  smaller  changes  in  the  heads  rate  for  the  same  analysis  after 
using  the  computer  than  before.   In  other  words,  the  operators  went 
through  a  learning  process  and  became  better  operators.   During  the 
period  when  the  computer  control  was  not  in  operation,  the  plant  was 
able  to  maintain  about  one  half  of  the  gains  obtained  with  the  computer. 

Other  advantages  were  also  apparent  after  the  installation  of 
the  computer  control  system  in  the  plant.   Some  of  these  were: 

1.    An  increased  reliability  in  instrumentations  due  to  the 
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TABLE  7 
Total  Unsaps  Concentration  in  Crude  Fatty  Acid 

Date  Time  *  Unsaps 

12-15-65  18:45  1.32 

22:30  1.46 

12-16-65  2:00  1.77 

6:15  1.72 

9:45  1.07 

14:15  1.83 

18:00  1.58 

23:15  1.58 

12-17-65  0:00  1.52 

1:45  1,42 

6:45  1.43 

8:15  1.43 

10:45  1.62 

15:30  1,43 

18:00  1.32 

12-18-65  4:15  1.42 

7:30  1.42 

11:30  1,18 

14:15  1.01 

19:00  1.21 

22:15  1.46 

12-19-65  2:45  1.42 

12:45  1.38 
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demand  of  the  computer  and  the  rapid  knowledge  of  instrument  mal- 
functions 

2.  Faster  knowledge  of  yields  and  throughput 

3.  The  ability  to  spot  differences  between  shifts  in  the 
operating  of  the  equipment 

4.  The  ability  to  provide  engineering  assistance  from  the 
main  plant  in  Jacksonville  in  case  of  processing  problems. 


CHAPTER  V 

CONCLUSIONS 

Several  specific  conclusions  may  be  made  on  the  basis  of  this 
work. 

1.  It  is  possible  to  write  a  theoretical  description  of  an 
actual  industrial  distillation  column  with  sufficient  accuracy  to  allow 
its  control  by  a  digital  computer. 

2.  The  technique  of  fitting  a  rigorous  mathematical  model  to 
actual  plant  data  can  be  used  to  obtain  estimates  of  physical  data  which 
might  otherwise  be  unknown.   In  other  words,  it  is  not  always  necessary 
to  have  all  physical  data  completely  defined  before  a  model  may  be 
used,  but  actual  operating  data  may  be  used  to  fill  the  gaps  in  our 
knovs^ ledge . 

3.  As  pointed  out  in  Chapter  III,  if  a  rigorous  model  is 
developed,  its  utility  in  situations  which  vary  widely  from  the  actual 
conditions  used  to  estimate  physical  data  will  be  greatly  enhanced. 
However,  if  plant  operating  data  had  been  collected  and  curve  fitting 
techniques  used  to  force  a  relationship  among  variables  thought  to  be 
important,  extrapolation  of  these  relationships  would  be  a  risky 
business . 

4.  The  computer  control  of  the  Port  St.  Joe  Crude  Tall  Oil 
plant  is  economically  feasible,  based  on  economic  data  prevailing  at 
the  time  of  the  test. 

5.  The  power  ratio  convergence  technique  has  proved  itself 
to  be  a  useful  method  in  the  solution  of  steady-state  distillation 
calculations  and  should  have  applications  in  other  types  of  trial  and 
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error  calculations. 

Because  of  the  author's  participation  in  a  pioneering  computer 
control  project  on  a  commercial  scale,  several  general  conclusions 
should  be  mentioned  which  derive  from  the  overall  effort  but  are  not 
primary  to  the  work  discussed  in  this  dissertation. 

1.  There  is  a  necessity  for  confidence  -  building  interaction 
between  man  and  machine.   It  is  important  for  computer  experts  to 
realize  that  some  sacrifice  of  machine  efficiency  may  be  necessary 

in  order  for  the  operating  people  to  have  a  clear  understanding  of  the 
reasons  behind  the  computer's  control  strategy. 

2.  In  any  project  involving  computer  control,  engineering 
analysis  and  computer  simulations  should  be  completed  before  computer 
hardware  and  additional  plant  instrumentation  are  specified. 

3.  The  process  industry  has  to  learn  to  rely  more  on  its  own 
people  and  less  on  the  computer  vendor  and  manufacturer.   This  is  not 
due  to  any  overt  failure  on  the  vendor's  part,  but  rather,  the  fact  that 
a  company's  problems  are  usually  unique  to  the  point  that  they  are  best 
defined  by  its  own  people. 

4.  The  reliability  of  computers  for  process  control  should  be 
one  of  the  first  considerations  in  specifying  the  computer  hardvvare. 
Control  applications  are  much  more  critical  than  traditional  data 
processing  to  machine  failure, and  downtime  due  to  computer  hardware 
must  be  held  as  low  as  possible. 

5.  The  installation  of  a  computer  on  an  industrial  process 
returns  many  benefits  that  are  difficult  to  justify  economically.   Such 
intangible  items  as  the  better  understanding  of  the  phenomena  taking 
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place  resulting  from  the  initial  study  of  the  process,  the  reliable 
and  up-to-date  information  available  once  the  computer  is  installed, 
and  the  ability  to  recognize  deviations  from  normal  operation  are 
items  which  cannot  be  given  a  specific  money  value,  but  are  neverthe- 
less important  results  of  such  a  project. 


LIST  OF  SYMBOLS 

Upper  Case  Letters 

A.  Antoine  Equation  Constant  (See  Equation  63) 
AN  Acid  Number 

B  Bottoms  Stream  Flow  Rate 

B.  Antoine  Equation  Constant  (See  Equation  63) 

C.  Antoine  Equation  Constant  (See  Equation  63) 
CFA  Sidestream  Flovv/  Rate 

D  Overheads  Flow  Rate 

E  Murphree  Liquid  Efficiency  (See  Equation  37) 

E  Murphree  Vapor  Efficiency  (See  Equation  14) 

F  Feed  Rate 

H  Total  Enthalpy  of  Organic  Vapor  Leaving  Stage  N 

H  Enthalpy  of  Steam  Leaving  Stage  N 

H.  Enthalpy  of  Vapor  of  Component  i 

K  Equilibrium  Constant  Defined  y  =  Kx 
n 

L  Liquid  Rate  Leaving  Stage  N 

N  Stage  Number 

i^  Rate  of  Heat  Gain  or  Loss 

R  Ideal  Gas  Constant 

RA  Rosin  Acid  Per  Cent 

RAN  RA  -V  1.855 

S  Steam  Rate  Entering  Reboiler 
o 

V  Vapor  Rate  Leaving  Stage  N 

XAN  Empirical  Acid  Number 
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Lower  Case  Letters 

a     Constant  in  Heat  Capacity  Equation  (See  Equation  65). 

b     Moles  Flow  Rate  of  Component  i  in  the  Bottom  Defined 
i 

b.  =  Bx, 
1      1 

b.    Constant  in  Heat  Capacity  Equation  (See  Equation  65) 

c    Total  Number  of  Organic  Components 

c     Heat  Capacity 
P 

d.    Moles  Flow  Rate  of  Component  i  in  CK-erhead  Defined 

d.  =  Dx. 
1     1 

e    Base  of  Natural  Logarithm 

f    Number  of  Feed  Stage 

g.    Gain  Factor  as  Defined  by  Equation  58 

h     Total  Enthalpy  of  Liquid  Leaving  Stage  n 
n 

h.     Enthalpy  of  Liquid  Component  i 

k.     Ratio  Raised  to  a  Power 

1 


in 


m 


Pure  Component  Vapor  Pressure  of  Component  i  at  Tempera- 
ture of  Stage  n 


t     Temperature 
X.     Concentration  of  Component  i  in  Liquid  Leaving  Stage  n 


Concentration  of  Liquid  vhich  \*)uld  be  in  Equilibrium 
with  Vapor  of  Concentration  y 

in 

y.     Concentration  of  Component  i  in  Vapor  Leaving  Stage  n 

y      Concentration  of  Water  in  Vapor  Leaving  Stage  n 

* 
y      Concentration  of  Vapor  Which  Would  be  in  Equilibrium 
with  Liquid  of  Concentration  x. 
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Greek  Letters 

a.  Relative  Volatility 

Q.  Factor  Defined  by  Equation  44 

X.  Latent  Heat  of  Vaporization 

■n  Total  Pressure  at  Stage  n 
n 

0  Arbitrary  Function 

-y.  Ratio  of  Composition  as  Defined  by  Equation  55 

Superscripts 

N  Kelers  to  Ne«/  iriai  value 

0  Refers  to  Previous  Trial  Value 

Subscripts 

B  Refers  to  Bottom  of  Tower 

C  Refers  to  Condenser 

c  Refers  to  Sidestream 

D  Refers  to  Overhead  of  Tower 

F  Refers  to  Feed  Stream 

f  Refers  to  Feed  Stage 

R  Refers  to  Reflux 

S  Refers  to  Reboller 

s  Refers  to  Steam 

T  Refers  to  Top  of  Tower 

w  Refers  to  Water  Vapor  in  Organic 
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APPENDICES 


APPENDIX  A 
EXAMINATION  OF  THE  DEGREES  OF  FREEDOM  FOR  A- 30  COLUMN 

An  actual  distillation  tower  is  a  determinate  system.   Under  a 
set  of  fixed  conditions  the  tower  produces  a  unique  set  of  stream 
compositions,  temperatures,  pressures,  and  flow  rates. 

The  purpose  of  any  mathematical  model  is  to  give  the  same 
solution  as  the  operating  unit  when  subiected  to  the  same  set  of  fixed 
conditions.   The  problem  of  vi/hich  variables  and  how  many  variables  may 
be  fixed  before  the  model  yields  a  unique  solution  has  been  considered 
by  many  authors  and  this  development  follows  closely  the  one  used  by 
Smith,  Design  of  E-iuilibrium  Stage  Processes,  Chapter  3. 

Essentially,  any  method  to  determine  the  degrees  of  freedom  of 

a  system  consists  of  counting  the  number  of  variables  and  the  number 

of  independent  equations  relating  these  variables.   The  difference 

between  these  two  quantities  is  the  number  of  variables  which  must  be 

specified  to  make  the  system  determinate  and  can  be  calculated  by  the 

simple  equation: 

N.^  =  N   -  N  (69) 

df    v    e 

In  a  complex  tower,  it  is  advantageous  to  consider  each  individual 

type  of  unit  separately  and  then  combine  these  elements  to  give  the 

system's  degrees  of  freedom. 

A.   Partial  Reboiler 
A  single  homogeneous  stream  will  contribute  C  +  2  variables  to 
any  element  or  unit  of  which  it  is  a  part.   If  the  stream  is  a  mixture 
of  C  components,  these  variables  are  the  C  -  1  compositions,  temperature. 
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pressure,  and  rate.   (Alternately  these  variables  may  be  considered 
the  C  +  1  degrees  of  freedom  as  predicted  by  the  thermodynamic  phase 
rule  plus  the  flow  rate.  ) 

For  a  partial  reboiler,  as  shown  in  Figure  34,  we  have  C  +  2 
variables  for  each  stream.   C  is  equal  to  the  number  of  organic 
components  plus  1  for  the  steam.   A  heat  stream  contributes  one  variable, 
its  rate.   Therefore  the  total  number  of  variables  is: 

N  =  4(C  +  2)  +  1  =  4C  +  9  (70) 

The  number  of  independent  equations  which  can  be  written  is 

C  material  balances,  1  energy  balance,  C-1  equations  relating  the 

liquid  and  vapor  compositions  leaving  the  reboiler,  and  2  equations 

relating  temperature  and  pressure  of  the  exit  streams;  the  total 

equations  are: 

N  =  2C  +  3  (71) 

e 

The  degrees  of  freedom  are: 

N^^  =  4C  +  9  -  (2C  +  3)  =  2C  +  6  (72) 

dt 

A  choice  of  variables  might  be  to  define  the  two  entering 

streams,  2(C  +  2),  the  temperature  of  the  exit  stream,  and  the  heat 

load,  q.   It  should  be  noted  that  in  specifying  the  steam  composition, 

the  organics  (C-1  of  the  components)  are  zero. 

B.   Total  Condenser 


The  condenser  system  on  A-30  distillation  column  can  be 
described  as  shown  in  Figure  35.   The  number  of  variables  is: 

N   =  3(C  +  2)  +  1  =  3C  +  7  (73) 
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Vapor  (Organic  and  Steam) 


Liquid 


^  Liquid 


Figure  34.   Schematic  of  a  Partial  Reboiler 


Steam 


Vapor 


Figure  35.   Schematic  of  a  Total  Condenser 
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The  number  of  independent  equations  ir.  : 

N  =  C  +  3  (74) 

The  number  of  independent  specifications  \vhich  must  be  made 
Is: 

N^^  =  2C  +  4  (75) 

These  variables  might  be  picked  by  defining  the  entering 
stream  (C  +  2)  and  the  temperature  and  pressure  of  the  condenser.   There 
are  C  variables  picked  to  show  the  lack  of  organic  in  the  steam  and 
the  lack  of  steam  in  the  liquid  organic. 

C.   Stream  Divider  (Reflux  and  Reboiler) 
A  stream  divider,  such  as  that  shown  in  Figure  36,  is  used  to 
split  one  stream  (L^ )  into  2  streams  (L„  and  L  ).   The  arrow  head  on  both 
ends  of  the  line  for  the  energy  stream  is  meant  to  show  that  there  may 
be  either  a  gain  or  a  loss  of  energy. 

The  number  of  variables  for  a  divider  is: 

N   =  3(C  +  2)  +  1  =  3C  +  7  (76) 

The  equations  which  may  be  ^^7ritten  are:   C  -  1  equation 
expressing  the  concentration  equalities  between  the  entering  stream 
and  each  of  the  leaving  streams  for  a  total  of  2C  -  2;  1  overall  material 
balance;  1  energy  balance;  and  2  equations  expressing  the  temperature 
and  pressure  equalities  between  the  leaving  streams.   Therefore: 

N  =  2C  +  2  (77) 

e 

The  number  of  independent  variables  which  must  be  specified  is: 

(78) 
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Figure  36.   Schematic  of  a 
Stream  Divider 
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Figure  37.   Schematic  of  a 
Equilibrium  Stage 
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These  variables  could  be  the  specifications  of  the  feed 
(C  +2),  the  heat  loss  or  gain,  the  pressure,  and  the  ratio  of  1.^11    . 
If  the  divided  stream  is  liquid  it  must  be  specified  that  the  composi- 
tion of  steam  is  zero,  which  only  leaves  C  +  1  additional  specification 
to  be  made  on  the  feed . 

D.   Equilibrium  Stage  (Fiaure  37) 
For  an  equilibrium  stage  there  are  4  entering  material  streams 
and  one  energy  stream,  giving  a  total  variable  count  of: 

N  =  4C  +  9  (79) 

As  independent  equations  we  can  \irrite  C  material  balances, 

1  energy  balance  equation,  2  equations  relating  the  temperature  and 

pressure  of  V   and  L  ,  one  equation  specifying  the  steam  composition 
m      m 

is  zero  in  the  liquid  phase  and  C  -  1  equilibrium  relation  between  the 

organic  components  in  L  and  V  .   The  number  of  equations  is: 
m      m 

N   =  2C  +  3  (80) 

This  leaves  the  degrees  of  freedom  at: 

N^^  =  2C  +  6  (81) 

If  the  stage  is  not  a  perfect  stage,  this  introduces  C  more 
variables,  the  composition  that  should  be  in  equilibrium  with  one  of 
the  exit  streams,  but  it  also  introduces  C  additional  equations,  the 
efficiency  equations. 

E.   Feed  Stage  (Figure  38) 
For  a  feed  stage,  there  are  5  entering  streams  and  one  energy 
stream  which  gives  a  total  variable  count: 


96 


^f+1 


->    ^ 


< >  Q 


f-1 


Figure  38.   Schematic  of  a  Feed  Stage 
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Figure  39.   Schematic  of  a  Sidestream  Stagt 
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N  =  5(C  +  2)  +  1  =  5C  +  11  (82) 

The  number  of  independent  equations  which  can  be  written  is 
the  same  as  for  an  equilibrium  stage: 

N  =  2C  +  3  (83) 

The  number  of  independent  variables  that  have  to  be  specified 
is: 

(84) 

These  specifications  are  generally  the  same  as  for  an  equilibrium 
stage  with  the  additional  C  +  2  variables  used  to  describe  the  feed. 

F.   Side  Stream  Stage  (Figure  39) 
The  number  of  variables  for  a  side  stream  stage  is  the  same  as 
for  a  feed  stage: 

N^  =  5C  +  11  (85) 

The  same  equations  can  be  written  expressing  the  total  equality 

between  S  and  either  V  or  L  : 

m     m 

N  =  3C  +  4  (86) 

e 

This  gives  the  number  of  independent  specifications  which  have 

to  be  made  as: 

(87)   ■ 

The  additional  variable  over  and  above  those  normally  made  for 
an  equilibrium  stage  could  be  the  side  stream  rate. 

G.   Combining  the  Individual  Elements 
The  degrees  of  freedom  for  the  total  system  may  be  calculated 
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by  combining  the  degrees  of  freedom  for  the  individual  elements,  taking 
into  account  any  decision  made  in  this  process  (such  as  number  of  like 
elements)  and  subtracting  any  restrictions  which  may  arise  from  this 
process . 

UTien  combining  K  equilibrium  stages  as  shown  in  Figure  40, 
the  number  of  variables  in  this  case  is  K  times  the  number  of  degrees 
of  freedom  from  each  stage  plus  1  to  reflect  the  choice  made  of  K 
elements: 

N  =  K(2C  +  6)  +  1  =  2KC  +  6K  +  1  (88) 

The  new  restrictions  which  arise  are  equations  which  express 
the  equality  between  streams  leaving  one  element  and  entering  adjacent 
elements.   There  are  2  (K-1)  such  interstreams ,  as  shown  in  Figure  40. 
This  gives  an  additional  number  of  equations: 

N   =  2(K  -  1)(C  +  2)  =  2KC  +  4K  -  2C  -  4  (89) 

e 

The  number  of  independent  variables  which  may  be  specified  is: 

N^^  =  2C  +  2K  +  5  (90) 

df 

A-30  distillation  column  is  shown  schematically  in  the  Figure 

41. 

Adding  up  the  degrees  of  freedom  for  each  element: 

Reboiler  2C  +  6 

Bottoms  divider  C  +  5 

Stages  below  sidestream  2C  +  5  +  2(S-1) 

Sidestream  2C  +  7 

Stage  below  feed  2C  +  5  +  2(M-(S)) 

Feed  stage  3C  +  8 

Stages  above  feed  2C  +  5  +  2(N-(M+1)) 

Condenser  2C  +  4 

Reflux  splitter  C  +  5 


Total  17C  +  46  +  2N 
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Figure  40.   Schematic  of 
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Figure  41.   Schematic  Representation 
of  A-30  Distillation  Column 
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There  are  14  interstreatns ,  which  leaves  as  the  number  of 
independent  variables: 

N   =  17C  +  2N  +  46  =  14C  -  28  =  3C  +  2N  +  18  (91) 

Due  to  the  nature  of  the  process  in  A-30  the  following  conditions 
must  be  met:   There  is  no  organic  material  in  the  steam  entering  the 
reboiler;  therefore  C  -  1  equation  may  be  written  specifying  zero 
composition  for  the  organic  components  in  this  stream;  the  assumption  is 
made  that  no  organic  material  is  lost  in  the  exit  stream  from  the  top 
of  the  column,  therefore  C  -  1  equation  reflecting  this  assumption  can 
be  written;  and  since  the  steam  is  assumed  to  be  present  only  in  the 
vapor  phase,  a  total  material  balance  on  steam  can  be  written.   This 
leaves  the  number  of  independent  variables  which  must  be  specified  before 
the  solution  to  the  equations  is  unique  as: 

(92) 

H.   Choice  of  Operating  Variables 

The  choice  of  these  variables  on  an  existing  column  is  usually 

made  as  follows: 

Pressure  at  each  stage  N 

Pressure  at  condenser  1 

Pressure  at  reboiler  1 

Pressure  at  2  splitters  2 

Heat  loss  at  each  stage  N 

Heat  loss  at  2  splitters  2 

Total  stages  1 

Feed  location  1 

Sidestream  location  1_^ 

2N  +  9 

The  above  considerations  leave  C  +  10  variables  to  be  fixed 

before  the  column  or  the  equations  are  determinate. 
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If  the  feed  is  specified  this  takes  C  +  2  of  these  variables  (with 
consideration  that  the  steam  composition  is  zero)  .   If  the  steam  input 
to  the  kettle  is  specified,  3  more  degrees  of  freedom  are  removed. 
This  leaves  5  variables  which  must  be  fixed.   These  variables  may  come 
from  the  following  list: 


Variable 

Reflux  Temperature 
Distillate  rate 
External  Reflux  Ratio 
Reflux  rate 
Condenser  heat  load 
Concentration  of  1  component 

in  one  stream 
Sidestream  rate 
Bottoms  rate 
Reboiler  Heat  Load 
Rate  to  Reboiler 


Degree  of  Freedom 


Some  of  the  above  variables  cannot  be  specified  indiscriminately. 
Some  are  restricted  to  narrow  ranges.   For  example,  the  condenser 
heat  load  must  result  in  a  condensate  uemperature  somewhere  between  the 
bubble  point  and  the  freezing  point  of  the  mixture.   Some  of  the  variables 
are  so  closely  related  that  it  is  impossible  to  specify  values  for  all 
of  them  independently.   For  example,  for  a  given  feed  rate,  the  distillate 
rate,  sidestream  rate,  and  bottoms  rate  must  give  an  overall  material 
balance  and  therefore  only  two  of  them  can  be  specified  independently. 


APPENDIX  B 
PHYSICAL  DATA  FOR  THE  TALL  OIL  FATTY  ACIDS  SYSTEM 

A.  Vapor  Pressure 

Vapor  pressure  data  for  the  various  tall  oil  fatty  acids  were 
collected  from  several  sources  (41,42,43')  and  available  industrial 
information  (44,45).   Figure  42  shows  data  for  oleic  acid  (C-18) 
which  Is  the  major  component  of  the  system.   Figure  43  shows  data  for 
palmitic  acid  (C-16).   The  data  for  stearic  acid  (C-20)  and  myristic 
acid  (C-14)  are  shovm  in  Figures  44  and  45.   An  available  Glidden  Co. 
report  (45)  gave  vapor  pressure  curves  for  three  unsaps  and  these  are 
shown  in  Figure  46. 

Using  the  vapor  pressure  data  as  discussed  above,  a  steady-state 
calculation  was  made  on  the  computer  and  the  resulting  column  profiles 
were  plotted  along  with  the  experimental  data  from  Run  040165.   The 
constants  in  these  curves  vjere  changed;  two  additional  unsaps  were 
postulated  to  correspond  to  current  analysis  information,  and  calcula- 
tions were  made  until  agreement  between  the  calculated  profiles  and  the 
experimental  points  were  within  what  was  estimated  to  be  the  experimental 
error.   The  final  constants  used  in  this  work  are  shown  in  Table  8. 

B.  Specific  Heats 

The  constants  used  for  the  specific  heat  curves  are  shown  in 
Table  9. 
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Figure  42.   Oleic  Acid  ^'apor  Pressur 
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Figure  43.   Palmitic  Acid  Vapor  Pressure 


100.0 


G     Reference  41 
^     Reference  42 


10.0 


1.0 


J \ \ \ \ \ ! \ L 


2.1 


2.2  2.3 

3 


2.4 


2.5 


t   +  230 


Figure   44.      Stearic   Acid   \'apor   P 
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Figure  45.   Myristic  Acid  Vapor  Pressure 
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TABLE 


Final  Antoine  Constants  for  Tall  Oil  Fatty  \cids 


logP. 


t  +  C. 


Component 

\ 

8.833 

A 

C-IA 

-3333.0 

230 

C-16 

-3247.0 

8.444 

230 

C-18 

-3802.3 

9.384 

230 

C-20 

-3590. 

8.8 

230 

L-1 

-3535 

9.42 

230 

L-2 

-3558.7 

9.319 

230, 

L-3 

-3755.0 

9.591 

230, 

L-4 

-3831,0 

9.610 

230, 

L-5 

-3852.0 

9.592 

230, 
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TABLE  9 


Specific 

Heat  Constants 

for  Tall 

Oil 

Fa 

.tty 

Acids 

Si 

=  a_ 

l^\' 

Component 

-122 

h 

Rej 

Eerence 

C-16 

0.9 

46 

C-18 

-  51 

0.9 

46 

C-20 

-198 

1.17 

46 

C-14 

-A- 

L-1 

226 

0.2667 

45 

L-2 

for 

t  < 

2120C 

L-3 

-  36 

1.5 

45 

L-4 

t  > 

212' 

^C 

L-5 

Mean  specific  heat.   Liquid  enthalpy  is  calculated  by  equation 


APPENDIX  C 
FORTRAN  LISTINGS  OF  CONVERGENCE  COMPUTER  PROGRAMS 


112 


vj   ~ 


UJ  _i  cci  "  o        — 
on:    -  o  -t    -  o 


a.  ly)  :^: 

ijj  n 
5:  I--  (- 
o  < 


O  IJ 1 

LU  IJJ 

(v  ^-  or 
a      a: 

OXo 

UJ  >-  q; 
~-   'I  UJ 

^-  ^  ^ 
3  ^  ^ 

o  ^ 

DC 

CO  U- 
^  O 
1/5 


CL   UJ 
O  D 


to  to  o 
^  X  o 

►  s:    ► 

'-I  O  f\i 

CL    2:    f- 
<t   LU   1/1 

oo  a  s: 
2:    ►  o 
-00 
UJ  5: 
o  o    ► 

<I     O     r-l 

t-  o  h- 


UJ   00 

UJ  c^  - 
U-  -  z 
i;  00  1—1 
^   ►  <t 

UJ  Q   O 

o    ». 

'/I  u.   < 

;c  o  f- 


o 

o    - 
-  > 

-•  X 


o  O 
o 


)  o 


X       ».  —    Q_ 

>  U-  □ 
•«  Ul  X  (_ 
U-  u_ 

X    -    -x 

0-.   — •     ^ 

■   ►  <j:  o  — 

Q.      I-      r-^     O 

O    l>0    —    — 4 

H-  Q.  cy-  — 

u.        01- 
X    -  q:  o 

Ql  Q'  CD 
-  -=I   UJ  U. 

^-  H-       X 

o  00    -  ». 

CO  > 

u_        oj  o 

X     ►  ^  <J 

o:  —  -^ 
•  <c  I-  _j 


-  o  — 

—  -^  _j 

O     -  X 

o-  o  - 

o_  —  o 

s:  or  -J- 

UJ  <  ,. 


C3    -  ::• 

-t  -   X 


•-  U.    00  UJ      ► 
-•-«         X  Q  '- 

m  >  ►  --a 
00  X  >  -  <^ 
^  2       >  2  - 

00-      o  > 

o  U  10     ►  W-, 
X   X  00 

UJ 


UJ     ^  — 

or  o  o 
a  -I  <f 


X  al 


■— I  <Ni  m  s}   in  ^c 


—  —  vO 

00m 


-•-    Q. 

a:  5: 
LU  o  - 
I-  (^)  f\j 
-)      t— I 

X  vO 
—I  (NJ 


O 
00 

~  <r  < 
UJ  -.  5; 

:?-  a:  cc 

O  U.  U- 


-  <  < 

[  or  a: 

■  o  o 

u.  u. 

in  vo 


LU    .-< 

-  <     + 


Q  (X 
—  <t  UJ 
U.  UJ  I— 
<-i  of   t-, 


u  O  O  o  o 


t\J  l-l      I-. 


113 


—        "        -o 


o 


(/5  ^.  ^  ^ 


o 

o 

fM 

< 

i^ 

II 

r>' 

— 

II 

o 

a: 

UL 

-.< 

o 

IJJ 

►-< 

;^ 

1—    M    (— 

1— 1  -^  — ~  w-  ,-, 

m  — .  r<^  — 

3      r) 

,.-,    ^^   q;;   ,_, 

O   —  O    ►-<   LLI 

o  *-•  o 

in  —  —  o  — 

r-t    ^    ^     —     3 

m   <\l    r-<    r^    o 

o       o  ;? 

LU   m  LU 

o  h-  1-  «:  1- 

o  1-  o  f-  '-' 

1—  rg  h- 

'-<  I/)    00   IJJ   (/I 

1-    LO  P-    00    1- 

t— 1        t— 1 

LU    UU   —  LU 

(u       lu  :2- 

q;  o  '"'^ 

o  1—   1—  U.  h- 

o  t-  O  f-  o 

o  q;  o  '"'^ 

O  3^  Q  ~ 


(>   O  fM  -J    in   O 


O  M  .-.  M  >-.  c'l  I-'  o  '-'  O 

fxi  f\i  f^J  fNj  f\i  r\j  rn  m  m  «^ 
o  o  o  o  o  o  o  o  o  o 


in 

UJ 

lU 



—  O   LU 

rg  o 

t- 

1- 

o 

O   vO    Q. 

o 

i-A 

iC 

;r  o  < 

ro 

in    ► 

+ 

+ 

1—1 

M  ^  1— 

•■ 

rg  —1 

< 

<    ► 

o 

— . 

,— 

o 

O  1^  1— 

rf> 

—    M 

^_| 

,_, 

in  3 

•. 

(n 

• 

— 

II 

o 

II   o  tx 

O 

\     ►"-! 

o 

o 

in 

■-<  1— 

o^ 

rn: 

1- 

f- 

-^ 

o 

--       3 

LU 

— 

LU  O 

00 

oo 

»-l 

r-l 

>-  —  o 

3 

o 

1-    vO 

UJ 

LU 

^^ 

^ 

z 

.-H    O 

H- 

I— 

Q 

o 

o  o  "JJ 

t-< 

1 

H-l 

>-i 

^ 

h- 

2  1-  H- 

H- 

^ 

_^ 

^^ 

-^ 

2 

— 

11-  o 

LL 

LL 

< 

o 

<  O  °i 

O 

U- 

-.  o 

"-^ 

•-' 

t5 

cr> 

O  O  3t 

o 

•-< 

a-  o 

_, 

rvj 

m 

<!- 

in  o  iri 

o 

in 

m  4 

<f 

■4- 

vf 

<^ 

<)■  in  in 

^r 

fM 

o  o 

O 

O 

o 

O 

O  o  O 

c- 

Q 

^ 

^ 

o  — 

■^ 

1— 1 

_ 

z 

in  ^-< 

'- 

O  1-1 

z 

„ 

_ 

S 

LU 

< 

r—l     

X 

<  — 

~. 

l-H 

I— 1 

1— 

LU 

7^ 

o 

—   CC 

-:r 

LL  z.\. 

in 

3 

~_    ,— 

— 

0 

2! 

X 

■}{■ 

-  O 

< 

in 

•^  O 

(>j 

LU 

3:    1— 1 

CL 

CO 

X 

(XI 

»• 

^{■ 

•"-I     1— 

u. 

r\i  — . 

^  I— 

r-A 

i;^: 

LU  — 

0 

>- 

-;;- 

D.    -5 

^ 

"  o 

o 

r-H  uj 

LU   >< 

•- 

X 

-i.  0 

f— 

3 

Q 

C3 

<I       - 

•- 

»-4 

>-H     < 

— 

-o 

(J  ;} 

CO 

-;!• 

X  X 

X 

>— 1 

— 

00    -0 

-0 

^^ 

—   LL 

1 

TO   < 

<  Q 

r-J 

,— 

\ 

00 

1 

2; 

Qi 

X  >:- 

.— 

O  1- 

h-   1 

.—1 

LU 

+  + 

+ 

— 

5; 

— - 

.-  LLI 

UJ 

o 

K-1 

to  — - 

».    — - 

0 

1—1 

3 

1—1 

"    Cl 

CL 

^ 

,— ( 

w- 

^  ;- 

^    >-H 

00     r-t 

»_4 

1-  Q 

O- 

~_ 

00 

— 

m  < 

<   0 

a: 

(_j  •—    ►■ 

LL 

CO     '•■ 

.»   _- 

00 

O  '-' 

0 

."S 

CL 

\ 

LL 

C\J  1- 

uu 

^  o  >-• 

X 

•— '  H. 

r- 1     1—1 

;^ 

moo 

LJ 

'JJ 

0 

■— 

X 

(NJ 

(_> 

•- 

(.-) 

o 

O 

O 

^  <  — 

;; 

,-4    — 

'—  >C 

_ 

•• 

~  7- 

Tu 

^ 

;^ 

1— 

1—1 

-  h- 

.<. 

h-    r-K 

;- 

II 

• 

• 

• 

-<  1-  >< 

LL 

—   X 

X    " 

—   X 

1— < 

:j  r) 

3 

•. 

X 

s: 

— 

II 

vO  r> 

•■ 

3 

•. 

o 

o 

o 

l>0 

I.L 

t-A 

--' 

OO  00 

00 

.— < 

3 

0 

rc^  CL 

■— 1 

CL     II 

,— 1 

_ 

II  2:  II 

II 

l_,    II 

Ln 

II    •— 

—       II 

LO 

X 

II 

II 

00 

X 

— ~ 

r-t    (— 

II 

1— 

fMIJJ 

II 

n 

II 

II 

II 

~^- 

tM 

II 

.5 

<M 

It    II 

II 

UJ 

^ 

1— < 

0 

t— < 

D  •-' 

3 

_ 

t-i  ^-. 

.— 

Q.  — 

r-4 

LU   — 

,— 1 

II 

3 

— 

— . 

II 

— ■ 

—  a 

0 

II   2 

2: 

q: 

a. 

1— 

Q 

O     —    H-< 

>-< 

o  — 

1— 1  1— 1 

z  •-* 

H-  Q 

CL 

z: 

ir> 

>-< 

1—1 

m 

00 

in 

1—. 

c<: 

_! 

a 

;-i 

■o 

r<1  X  — 

1 — 

o 

.^    w. 

^<  — 

0 

-~ 

Q  '-' 

0 

v-l 

m 

-- 

— 

-^ 

UJ 

0   LU 

f<^  LU  m 

ei- 

3 

O 

i-- 

CO 

oo 

'-*        s. 

D. 

X  (X 

1-  ql  3: 

:« 

\^ 

>— 1 

OD  10 

1— 

1— 

3 

CL 

»-l 

CL 

1-    1- 

1— < 

1-  fv 

oe 

(— 

o 

"• 

jr_ 

s: 

—  LU 

O 

-'  o 

O  UJ 

—  LU 

— 

s:  s: 

s: 

z: 

ID 

0 

— 

7- 

1— < 

1— • 

oc 

i.U 

o 

< 

5 

Z) 

3 

O  LL   2 

1- 

LL  H- 

s 

1-3: 

IL   2: 

0 

0 

Z)  =) 

o 

0  0 

Z 

H- 

0 

0 

0   ^- 

0 

a-.  0 

-)'Oo^ 

o 

li. 

(/) 

I/) 

1/1 

O  -<   X 

X 

'-•  X 

X  X 

•-H     X 

0 

X 

00  00 

00 

0  0 

X 

X 

X 

0 

0  3 

0 

3  0 

o  o 

o 

o 

^ 

eg 

ro 

vl 

O    -•   LO 

O  (^  00 

cr 

in  vO 

r-  CO 

cr 

Ln 

vo  r- 

00 

0 

_, 

f\i  m  -4- 

in 

vO 

t-- 

CO  in 

f.  -t 

cr 

o^ 

(T- 

cr 

a> 

o  o  o 
— <  —1  •-< 

o  o  o  o 

r-l 

f\i  (\j 

rj 

r<^  ro 

m 

m 

ro 

ro 

r^ 

ro  rvi 
,-1  nj 

115 


►-■  o  —  o  3  I-  Q 

—  o  a.  u  a  UJ  ^ 

><")>-<  -7  o  a  UJ 

m  c  CM  m  in  o 

f^^>  -;)    vf  -i  -t  ir\ 

r-i  oj  oj  r^j  f\J  rvj 


116 


I-  a-  I  o 

_j  —  »-   ► 

X  a  t-  o 

-o  D    ► 

'-  H-  CO  00 

a  U 1  o 

—  X  X  < 


CT'  CO 
—  LL 


Ch  _ 

LL   h- 

~cr 

X   1- 

»• 

U-  — 

-  O 

X   D. 

U-   rO 

cr 

c^ 

►  o 

X      ». 

—  1- 

-  LL 

fNJ 

3: 

cr  > 

tx    ► 

1— 

UJ 

—  X 

OO 

00 

Q. 

H-  U.I 

(— 

X 

o   ► 

X  IJJ  1/1 

1—1 

1— 

►  IL 

•. 

X  — 

O     ►     - 

^ 

-^ 

►  c^ 

>—  a.  :^: 

—~ 

c> 

(/I  o  < 

a- 

r\i  — 

0>  Q 

X  1-  m 

h-  or 

—  2: 

»•>  (- 

,-< 

(.^  H- 

o.-< 

I-   X  00 

1- 

H-   (.) 

H-.  <T 

o    ► 

0-) 

1-1   < 

•-O  O    CO   (V      ^ 

h- 

LL 

X     ► 

X  UU  Q. 

i-j 

». 

,-, 

..  -^ 

j:*-  a 

•- 

r— 1         •- 

— 

—  0^ 

■ — 

H-   — 

CO  " 

Ov  — 

•— •    •■  u_ 

O- 

OO  0^ 

•  .~ 

—  o 

t/5  O.   X 

t—  — 

o  o 

(—  ■-• 

x  o 

o 

l-H     Ql 

r\j  r-l 

OOO 

.^h:;  ► 

1— 

o 

LU  -- 

D1  _) 

— <  I/)  1— 

00 

•-  q: 

Lr\  m 

X  X 

c-^  X  o 

H- 

o  <^ 

-^  ~- 

X  ^n 

— 1 

t— • 

p;  lij 

o  o 

OJ  O  IX 

I/O 

2^ 

►-  ^  X 

1— 

2! 

\-  e 

l-H    t-H 

o 

OO 

00 

a 

•-•  o 

»_    — 

1  - 

O 

>-< 

•-< 

lyl 

e  X 

.  00 

21  00 

1-    (- 

^ 

o 

;z 

O  ^ 

<     < 

Ul 

:r 

LU 

---   IJU 

?:  s: 

x 

:£ 

^ 

s:  z 

or  o: 

t— < 

o 

•— • 

o  >- 

o  o 

o 

o 

O  O  Q  lU  LL 

o^  o  q: 


•4-  vo  r-  00  o  •-'  Qi 

II  •    O  lU 

II    II    II    no    •  t— 

r-l    —    ^     —     —      II         II 

—  c\j  m  vj-  ir\  II    II 

—  —  —  --  CL    <  o 

Q.Q-Cl.Q-CL<tt—   nfO. 
•-I<t<t<<00       ILUO" 
l/IOOOOOO    002    (Uf-o 
-}-J-3-3-)r)0'-"       (J^ 


•-.  m  +  r- 


O  or   — 


OO 

LU  LU      •> 

00  3  -"I 

LU  H-i  I- 
00  I--   2 


1- 

p 

1 

OO 

•^ 

1— 

t- 

— 

'-' 

•— 

*~' 

M 

II 

a 

1—1 

1—1 

IX 

— 

-— 

LU 

(N) 

"-I 

— 

O^  O  (M  -J    in  o 
l-H  ^  --<  -4  r\) 


LA  vO  f^ 
f\l  r\J  tM 
O  O  O 


117 


v)   in 
—I  o 


4-  in 
o  -<f 
-•  o 


+ 
+ 

C\J 

iir   r\l 


0 

►_,  , 

.-1 

_ 

< 

„ 

(X 

0  in 

1— 

X 

in        '- 

^  !z! 

^ 

Q 

0 
m 

0 

1;^ 

0 

~  r\j 

0 

< 

0 

(\J                >-4 

►. 

1—     r-l 

>— < 

1— 

< 

K 

LL 

s 

— 

in 

1/1    1- 

'— 

<t 

♦-in 

— 

0 

U- 

-—  CO 

\ 

00          LU 

1— < 

ni 

t—     i/l 

1—1 

0 

0  ro 

t-> 

CO 

Jt 

>->  0 

^- 

■— 

'^ 

•— t          PZ 

— 

0^     1-4     h- 

— 

m    - 

— 

X 

— 

—  .— 1 

>— 1 

(— 1 

1— 1 

■-  1  —   X 

m 

•- 

l-l 

0 

>f 

0  '- 

a: 

(—4 

U-     •■ 

— 

— 

— 

-.-.  If 

1— 

in 

-     +       + 

-o^  ^- 

0 

,— ) 

C^     1 

_- 

><  00 

n. 

(X 

LL 

00  —  — 

I/) 

f\i 

^-1 

>— 1 

— 1 

0     »•  Q- 

rp  cc 

\— 

i;-  0 

0 

0 

X 

H- 

.    .~ 

<t 

* 

o>  in  h- 

oc. 

+ 

►  ^- 

CD 

IL     r-l 

1- 

h- 

■;:- 

r-l    p-|     • 

^ 

0 

•-• 

f^ 

11    y-i  f-> 

'.D 

— 1 

(■0  =) 
-^  "  0 

j-ai 

0 

0 

0 

Q. 

r-,    !_. 

>< 

"    ~ 

X 

^ 

LL 

rn  Q 
—  X   •-• 

II  ir-i 

II 

•0 

II 

\ 

*-'  0  n 

11 

0 

II 

o>  — <  — - 

II 

• 

0 

II  a. 

II      t-> 

II 

II 

'^         oo 

m 

rci 

fy: 

h-  1— 

m 

1    ^. 

II 

0 

0 

0 

0 

0 

II 

~_ 

in 

II 

•-<    II    X 

—  0 

-^ 

0 

-^ 

UJ 

ai 

00  00 

, — 

0 

— 

UJ 

1        2: 

--^ 

-.  H 

— ~    1— 

— 

(Nl 

-^ 

~- 

1— I  —1 

•— < 

^-1 

t-i 

n 

h- 

0   (-  (- 

1-4 

>--i 

3 

oi  >— 

1— <  ►-< 

II 

II 

II 

X 

-— 

•-"   IX 

t— 1 

r-i 

t— 1 

'— 

jZ    —     \\ 

— 

— 

— 

— 

1— 1 

m  1— <  1— 1 

s_ 

^- 

:y<  UJ       , 

~- 

II 

0 

1-1 

■—  0 

— 

— 

t— 1 

LU  •-< 

a  n 

0 

0 

r-1 

>-j 

— 

0  —  — 

Q  0  Q 

1— < 

_-  h-   ' 

-f  a: 

Q- 

H 

0 

ro 

— 

h- 1 

1- 

0 

1— 

— 

^  — ^  ■— 

H-l- 

I-- 

H^ 

h- 

.-< 

1- 

-~ 

1- 

•— '      1 

O-   h- 

0 

0 

1— < 

Q. 

r-< 

3: 

D-   X 

Q. 

h- 

CL 

7^ 

X   .^    ►-4 

UO 

C/0 

oo 

2: 

(—1 

t-i 

p» 

—      , 

0 

t 

CO 

00 

0 

_ 

UJ 

0  — 

0 

0 

III 

—  UJ   — 

H-  0 

t— 

0 

^- 

0 

LL 

0  IJ-  li- 

< 

s 

<t 

□ 

IL  U.    •5t 

0  < 

^ 

;c 

0 

0  LL 

z:  t-u- 

h- 

s 

1- 

zz 

LL  e   0 

.-0 

•-* 

0 

*-" 

0 

0  --^ 

0 

0  U 

•-■»-<   0  0  IL 

00 

00 

00 

Q-. 

X 

X  •-" 

X 

X 

X 

»-i  X    X 

TO  O' 

0 

^ 

vh 

in 

c^ 

0   --•  (NJ 

m 

<f  in 

0 

ir\  0  in 

0  -1 

(M 

(^ 

-d- 

in 

0  --• 

in 

0  r- 

00 

in 

vO  r-  00  in 

(M(M 

r<^ 

m  m 

fO 

x^  -J-  vf 

<f 

<f 

1- 

in 

f\)  f<^  ro 

CT-  a- 

CP 

a- 

a- 

(P 

0  0 

0 

0  0 

0 

r-l 

--< 

-i^  rvj 

00 

r-l   r-1 

0 

0 

0 

0 

0000 

0 

0  0 

0 

-<r-( 

—1    r-l 

'"' 

•^ 

'"' 

— 4   i-l   — t 

lU 


r--        Q       — 


'Ir.R 

OO 

_ 

o 

>-i    .. 

o 

^, 

j: 

2 

2 

o< 

LJ 

ir\ 

•— < 

"-> 

1—1 

~j  — 

<t 

o  — 

00 

:3 

_) 

O- 

■— ' 

" 

XXX 

;:; 

h- 

00 

:5: 

-  vt 

^ 

X   -J 

< 

ir 

■>;- 

(X 

„ 

o 

■"* 

o 

+  +  + 

;5 

I/O 

— 

tn 

+  (J 
1— 

^ 

^ 

o 

lyi 

^ 

LL 

1^ 

1—  Q  (X 

••   LU 

1- 

(M 

•. 

a'.  >-H 

o; 

3 

LJ 

a: 

■it 

CO  '^  1— 

>< 

,_, 

f-4        »• 

—J 

,-, 

:k  .5 

LU 

1 

1 

\ 

(— 

LU 

^ 

— 

v^ 

„ 

—  in 

1 

o 

^-* 

00  ly) 

l~ 

n' 

Q 

>— < 

(- 

•-^ 

li. 

II 

It 

(_) 

(NJ 

•  II 

_ 

-- 

,—, 

q: 

fV 

;^ 

1 

:?: 

<— 1 

• 

U_ 

o 

-;;- 

Coirt  J) 

II 

X 

rvi 

D. 

o 

Q- 

-:t    LU 

3 

3 

OO 

1/1 

o 

Vr 

^ 

rO  IX 

LU 

,^ 

LU 

O 

o 

< 

OO 

LU 

•. 

LU 

00 

00 

;^- 

1 

UJ 

lU 

LU 

»• 

IJL) 

vO 

'— ' 

•   i!- 

II    11     II 

ID 

,^ 

,_, 

II 

3 

O  ~ 

o 

II 

00 

II  :^. 

3 

.—1 

3 

p^ 

z: 

< 

LL 

"^ 

00 

O 

■— t 

=:> 

1 

"— ' 

'—' 

O  rn 

2 

ir\  H-t 

,^ 

^ 

ro 

—) 

u.\ 

- 

z. 

< 

< 

^!• 

■z. 

2! 

z 

2r 

• 

1      • 

1—  Q  Q. 

m  — 

1— 

_ 

l_l 

O 

- 

or  s^ 

a:  00 

1— 

Q 

II 

*-* 

LU 

•— ' 

1— 

•—I 

Q 

o 

o 

Q  O 

CD  (/)  (— 

1— 

r-i  3: 

Q. 

,_, 

t— 

_'  •"" 

.^ 

s  — < 

II 

3:  -^ 

t- 

Z  H- 

^^ 

— 

(— 

OO 

1— 

71 

h- 

1— 

LU 

o 

. 

z: 

"i 

00 

00 

^ 

2: 

II 

< 

5r 

— 

z: 

•-» 

s: 

■ — 

II 

—    II 

>Z  ~.   ^ 

o  ("1  :s: 

1— 

o 

o 

^n 

LL 

^:  o 

-> 

Z:  LL 

O 

o: 

o 

UL 

U- 

IL 

o  u. 

o 

CC 

o 

LL 

Q 

LL 

(/loo  in 

(J 

5  X 

•X 

X 

o  o  o 

<  Q 

-5  <  i-H 

o 

Q. 

o 

•-* 

>-< 

o 

O 

LJ 

•-• 

O  Q- 

LJ 

•"■ 

a  u) 

<-'   O 

■or-  CO 

o 

.-1  f\j 

m  <»• 

U-N 

_, 

o  -< 

rvl 

rn 

r- 

o 

-4- 

ir\ 

o 

r-J 

(\J 

in 

4- 

rvj 

in  o 

rj  rv|  f\) 

m 

r*i  m 

r<^ 

r^ 

m 

J 

4  -J 

-t 

~t 

o- 

Ch 

<^ 

4- 

in 

in 

ITi 

o^ 

cr 

in 

in  in 

,— ) 

<-<  1— ( 

1— 1 

r-l 

1— t 

— 1 

•-'  r—l 

1—1 

o^ 

cr 

■"• 

<— < 

■"* 

'-• 

■"■ 

cr 

a- 

■"' 

.—1  .—1 

119 


^  ro 

t\J           „   ^_    ^  „  ^ 

Q,       >-i  —  rg  -I  I- 

o^  —  F-  a.  Q. 

—        3:  Q.  —  ni  h- 

fM     -UU   Q  sO  (—  3 

--t  ^  I-  rn  I/O  □ 

I  XXI 


U.I  ^ 


,  >f\  f-<  ^1  <-,  I  , 

r<-|  ^  —  ^  • 

— -  CM  t—  Q_  )  , 

O  O  '  t  I 

LL    n  'J3  I—  U-  <t  '^ 

•-  O  X  X  M  o  o 


iri  Lfi  a  >r\ 
ir\  r\j  m  ro 
r-<  rvJ  rg  rg 


m  o 


APPENDIX  d' 

GAS  CHROMATOGRAPHIC  ANALYSIS  OF  TALL  OIL  FATTY  ACIDS 

The  following  are  the  specifications  for  the  VPC  analysis  of 
Tall  Oil  I^atty  Acids. 
INSTRUI'lENTS: 

1.  Aerograph  (Manual  Temper^iture  Programmer  Gas  Chroma tograph) 
Wilkins  Instrument  and  Research,  Inc.,  Model  A-90-P. 

(1)  Thermal  Conductivity  Detector 

(2)  Isothermal  Operation 

2.  Recorder  (E.  H.  Sargent  and  Co.)  Model  SR 

(1)  Range  1 . 0  mv . 

(2)  Sensitivity  -  maximum  obtainable  while  maintaining 
satisfactory  stability. 

3.  Integrator  (Disc.  Chart  Integrator) 

Disc.  Instrument,  Inc.  Model  K4-1,  for  Sargent  SR  and  MR 
Recorders . 

4.  Chart  (E.  H,  Sargent  and  Co.)  Stock  No.  S-72168 

(1)  Scale  -5  to  100 

(2)  Zero  Right,  Integrator  scale 

(3)  Range  Plug  S-72181-30 

OPERATING  CONDITIONS: 

Carrier  Gas  -  Helium  100/140  ml/min. 

Temperatures: 


The  work  discussed  here  was  done  by  Mr .  G.  Watkins  of  the  Port  St.  Joe 
Lab.,  The  Glidden  Company. 
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Detector  280  -  290  C 


Injector  290°C 


Collector  220  -  230  C  (approx.) 

Filament  Current  -  200  ma. 

Attenuator 

S96C  (Methyl  Esters)  -  Initial  Setting  1,  Att.  C18  peak  to  8, 

return  to  1  for  remaining  peaks. 

A30  Feed  (Methyl  Esters)  -  Initial  Setting  2,  Att.  C18  peak 

to  8,  thence  to  1  for  remaining  peaks. 

S96C  Unsaps .  (Extracted)  -  Initial  Setting  1  (No  Attenuation) 

A30  Feed  Unsaps.  (Extracted)  -  Initial  Setting  1  (No  Attenuation) 

Chart  Speed  -  18"/hr. 

Sample  Size  - 

S96C  (Methyl  Esters)  8.0  Microliters 

A30  Feed  (Methyl  Esters)  7.5  Microliters 
A30  Feed  Unsaps.  (Extracted)  4.0  Microliters 
S96C  Unsaps.  (Extracted)        3.0  Microliters 

COLUMN; 

Length   -    4   Ft .    (copper   or   S.S.) 

ID   -    4mm 

Liquid  Phase  -  Apiezon  "L"  (13.77=  loading) 

Support  -  Chromosorb  "W"  (80  -  100  mesh  -  acid  washed) 

PROCEDURE  FOR  PREPARATION  OF  APIEZON  "L"  COLUMN: 

Weigh  1.4  grams  of  Apiezon  "L"  in  50  ml.  beaker.   Dissolve  in  15  ml.  of 
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chloroform  stock  solution.   Pour  solution  over  33  ml.  (packed)  of 
chromosorb  W  (.-jcid  washed).   Rotate  the  damp  mass  (don't  grind  or 
stir).   Remove  excess  chloroform  with  water  aspirator  at  room  temperature 
until  dry  (to  smell).   Pack  in  4  ft.  length  of  quarter  inch  copper 
tubing.   Plug  ends  with  glass  wool.   Flow  helium  through  column  while 
heating  for  first  time,  but  do  not  connect  to  detector  until  after 
heating  at  250   for  30  minutes. 

PROCEDURE  FOR  PREPAR/^TION  OF  CHLOROFORM  STOCK  SOLUTION: 

Dissolve  1  gram  Alkaterge  T,  1  gram  Tergitol  N.  P.  35,  and  1  gram  Span 

80  in  300  ml.  of  Chloroform. 

PROCEDURE  FOR  PREPARATION  OF  SAMPLES  FOR  G.C.  ANALYSIS; 
Fatty  Acids  -  Measure  approximately  10  -  12  drops  of  fatty  acid  sample 
into  a  50  ml.  beaker.   Add  5  to  10  ml.  of  Diazomethane  and  evaporate  on 
steam  bath  under  safety  fume  hood.   ASTM  D-1983-64T  BF  Methanol  method 
for  esterif ication  may  also  be  used  as  an  alternate. 

Unsaponif iables  -  This  method  is  applicable  to  all  products  derived 
from  Crude  Tall  Oil.   The  purpose  of  the  method  is  to  obtain  a 
representative  sample  of  unsaps .  for  gas  chromatography  analysis  as 
rapidly  as  possible,  therefore  all  measurements  are  performed  volumetri- 
cally . 

Procedure:   Pour  20  ml.  of  sample  into  250  ml.  erlenmeyer  flask.   Add 
40  ml.  of  alcoholic  K.O.H.  and  boiling  stone.   Reflux  on  hot  plate  with 
condenser  attached.   At  the  end  of  30  mins .  cool  flask  with  cold  water 
from  tap. 


123 


Transfer  to  500  ml.  separatory  funnel  using  100  ml.  distilled  water  to 
facilitate  transfer.   Extract  with  125  ml.  ether.   Add  ether  extract 
to  second  500  ml.  separatory  funnel.   Extract  again  with  30  ml.  ether. 
Add  second  ether  extract  to  second  funnel.   Wash  combined  ether  extract! 
with  15  ml.  distilled  water.   Evaporate  ether  extract  on  steam  bath  in 
100  ml.  beaker  containing  boiling  stone.   Ehry  unsaps .  in  oven  for 
5-15  mins .  as  time  permits. 

GENERAL  PROCEDURE  FOR  CALCULATION  OF  ALL  PEAK  AREAS: 

The  area  per  cent  of  each  peak  is  calculated  by  dividing  the  number  of 
integrator  strokes  under  each  peak  by  the  total  number  of  integrator 
strokes  for  the  entire  chromatogram. 

Area  per  cent  of  each  peak  is  converted  to  relative  weight  per  cent  by 
multiplying  each  area  per  cent  by  its  assigned  molecular  weight.   The 
total  molecular  weight  is  determined  by  the  summation.   Each  molecular 
contribution  is  divided  by  the  total  molecular  weight  to  obtain  values 
quite  close  to  true  weight  per  cent. 

PROCEDURE  FOR  CALCULATION  OF  G .  C.  ANALYSIS  FOR  COMPUTER: 
A- 30  TO^^fER  FATTY  ACIDS 
SYLFAT  96C  See  Figure  47. 

M.  W.  ASSIGNMENTS 
Peak  Area  #1               C14  195 

Peak  Area  #2  C14  234 

Peak  Area  0^3  C16  256 

Peak  Area  #4  C16  266 
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M.  W.  ASSIGNMENTS 

Peak  Area  #5 

CIS 

284 

Peak  Area  #6 

C20 

297 

Peak  Area  #7 

Unsaps 

315 

A-30  FEED  See 

Figure 

48. 

M.  W.  ASSIGNMENTS 

Peak  Area  #1 

C14 

195 

Peak  Area  #2 

C14 

234 

Peak  Area  #3 

C16 

256 

Peak  Area  #4 

C16 

266 

Peak  Area  #5 

C18 

284 

Peak  Area  #6 

C20 

297 

Peak  Area  #7 

L-3 

312 

Peak  Area  #8 

L-4 

318 

Peak  Area  #9 

L-5 

332 

n.  C.  ANALYSIS 

EXTRACTED  UNSAPS   CALCULATION: 

S-96C  AND  A-30F   See  Figures  49  and  50. 
Peak  Area  #1  L_  i 


Peak  Area  #2 
Peak  Area  43 
Peak  Area  #4 
Peak  Area  #5 


L-2 
L-3 
L-4 
L-5 


239 
290 
312 
318 
332 


Peak  Areas  are  converted  to  weight  per  cent.   This  weight  per  cent  is 
divided  by  the  calculated  unsaps    per  cent  to  obtain  actual  weight 
per  cent  of  unsaps   distribution  in  fatty  acid  sample. 

Fatty  acid  samples  are  corrected  for  unsaps.   Content  as  follows: 
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Figure  47.   Chromatogram  of  Crude  Fatty  Acid 
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Figure  48.   Chromatogram  of  Feed  to  A-30 
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Figure  49.   Chroraatogram  of  Unsaps  in  Feed  to  A-30 
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Figure  50.   Chromatogram  of  Unsaps  In  Crude  Fatty  Acid 
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C14  *  Less  1/2  L-1 

C16  Less  1/2  L-1 

C18  Less  1/2  L-2 

C20  Less  1/2  L-2 

CALCULATION  OF  UNSAPS  PER  CENT  IN  S-96C 


[  '     '  +  RA/o  =  Unsaps,  ( />) 


1°«     .   1.995  y 


RAN  =  RA7o  X  1.855 


CALCULATION  OF  UNSAPS  PER  CENT  IN  FEED 


XAN     ' 


RAN  =  RA7o  X  1.855 

XAN  derived  from  Figure  51 
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